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Hard ball systems are completely 
hyperbolic 

By Nandor Simanyi and Domokos Szasz* 
Abstract 

We consider the system of N (> 2) elastically colliding hard balls with 
masses mi, ... , rajy, radius r, moving uniformly in the flat torus Y V L = M u /L ■ 
W , v > 2. It is proved here that the relevant Lyapunov exponents of the 
flow do not vanish for almost every (N + l)-tuple (mi, . . . , m^v; L) of the outer 
geometric parameters. 

1. Introduction 

The proper mathematical formulation of Ludwig Boltzmann's ergodic hy- 
pothesis, which has incited so much interest and discussion in the last hundred 
years, is still not clear. For systems of elastic hard balls on a torus, how- 
ever, Yakov Sinai, in 1963, [Sin(1963)] gave a stronger, and at the same time 
mathematically rigorous, version of Boltzmann's hypothesis: The system of an 
arbitrarily fixed number N of identical elastic hard balls moving in the v-torus 
Y v = W /IS 1 (v > 2) is ergodic — of course, on the submanifold of the phase 
space specified by the trivial conservation laws. Boltzmann used his ergodic 
hypothesis when laying down the foundations of statistical physics, and its 
various forms are still intensively used in modern statistical physics. The im- 
portance of Sinai's hypothesis for the theory of dynamical systems is stressed 
by the fact that the interaction of elastic hard balls defines the only physical 
system of an arbitrary number of particles in arbitrary dimension whose dy- 
namical behaviour has been so far at least guessed — except for the completely 
integrable system of harmonic oscillators. (As to the history of Boltzmann's 
hypothesis, see the recent work [Sz(1996)].) 

Sinai's hypothesis was partially based on the physical arguments of 
Krylov's 1942 thesis (cf. [K(1979)] and its afterword written by Ya. G. Sinai, 
[Sin(1979)]), where Krylov discovered that hard ball collisions provide effects 
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analogous to the hyperbolic behaviour of geodesic flows on compact manifolds 
of constant negative curvature, exploited so beautifully in the works of Hed- 
lund, [He(1939)] and Hopf, [Ho(1939)]. 

The aim of the present paper is to establish that hard ball systems are, 
indeed, fully hyperbolic, i.e. all relevant Lyapunov exponents of these systems 
are nonzero almost everywhere. Our claim holds for typical (N + l)-tuples 
(mi, . . . , mjv ; L) € IR+ +1 of the outer geometric data of the system; that is - 
in contrast to earlier results — we do not require that particles have identical 
masses, though, on the other hand, we have to exclude a countable union of 
proper submanifolds of (N + l)-tuples (mi, . . . , mjy; L) (which set is, in fact, 
very likely to be empty). 

Full hyperbolicity combined with Katok-Strelcyn theory (see [K-S(1986)]) 
immediately provides that the ergodic components of these systems are of pos- 
itive measure. Consequently, there are at most countably many of them, and, 
moreover, on each of them the system is K-mixing. Our methods so far do not 
give the expected global ergodicity of the systems considered (an additional 
hypothesis to provide that is formulated in Section 6). 

The equality of the radii of the balls is not essential, but, for simplicity, 
it will be assumed throughout. For certain values of the radii — just think 
of the case when they are large — the phase space of our system decomposes 
into a finite union of different connected components, and these connected 
parts certainly belong to different ergodic components. Now, according to the 
wisdom of the ergodic hypothesis, these connected components are expected to 
be just the ergodic components of the system, and on each of them the system 
should also possess the Kolmogorov-mixing property. 

Let us first specify the model and formulate our result. 

Assume that, in general, a system of N(> 2) balls, identified as 1, 2, . . . , N, 
of masses mi, . . . ,mN and of radius r > are given in T U L = W /L ■ Z", the 
i/-dimensional cubic torus with sides L (y > 2). Denote the phase point of the 
ith ball by (qi,Vi) £ T V L xR". A priori, the configuration space Q of the N 

balls is a subset of f^' u : from T^'" we cut out cylindric scatterers: 

(1.1) C iJ = {Q = (q 1 ,...,q N )eT^ :\\q i -q j \\<2r}, 

1 < i < j < N, or in other words Q := T% u \ \Ji<> t<j <N^i,j- The energy 
H = i niivf and the total momentum P = J2i m iVi are first integrals 
of the motion. Thus, without loss of generality, we can assume that H = \, 
P = 0. (If P 7^ 0, then the system has an additional conditionally periodic or 
periodic motion.) Now, for these values of H and P, we define our dynamical 
system. 
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Remark. It is clear that actual values of r and L are not used, but merely 
their ratio L/r, is relevant for our model. Therefore, throughout this paper we 
fix the value of r > 0, and will only consider (later in Sections 3-4) the size L 
of the torus as a variable. 

In earlier works (cf. [S-Sz(1995)] and the references therein), where the 
masses were identical, one could and — to obtain ergodicity — had to fix a 
center of mass. For different masses — as observed in [S-W(1989)] - this is 
generically not possible and we shall follow a different approach. 

The equivalence relation ^ over Q, defined by Q ~^ Q* if and only if 
there exists an a G T V L such that for every i € [1, N], q* = qx + a allows us 
to introduce Q := Q/^ (the equivalence relation means, in other words, that 
for the internal coordinates qi — qj = q* — q* holds for every i,j £ [l,N]). 
The set Q, a compact, flat Riemannian manifold with boundary will actually 
be the configuration space of our system, whereas its phase space will be the 
ellipsoid bundle M := Q x £, where £ denotes the ellipsoid X)£Li m i(<^Zi) 2 
= 1, YliLi m idqi = 0. Clearly, d := dimQ = Nv — v, and dim£ = d — 1. 
The well-known Liouville measure fi is invariant with respect to the evolution 
S* := {S t : t G R} of our dynamical system defined by elastic collisions of the 
balls of masses mi, . . . , mjy and their uniform free motion. Here we have two 
remarks: 

(i) The collision laws for a pair of balls with different masses are well-known 
from mechanics and will also be reproduced in the equation (3.9); 

(ii) The dynamics can, indeed, be defined for \x — a.e. phase point; see the 
corresponding references in Section 2. 

The dynamical system (M, S*, n)m,L is called the standard billiard ball system 
with the outer geometric parameters (fh;L) £ 

Denote by i?o := Rq(N,v,L) the interval of those values of r > 0, for 
which the interior Int M = Int Q x £ of the phase space of the standard billiard 
ball flow is connected. From we should also exclude an at most countable 
number of values of r where the nondegeneracy condition of [B-F-K(1998)] (see 
the Definition on page 697) fails. The resulting set will be denoted by R$. 

The basic result of our paper is the following: 

Main Theorem. For N > 2,u > 2 and r G R , none of the relevant 
Lyapunov exponents of the standard billiard ball system (M, {S*}, fi)^ l van- 
ishes — apart from a countable union of proper analytic submanifolds of the 
outer geometric parameters (rh;L) £ 

An interesting consequence of this theorem is the following 
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Corollary. For the good set of typical outer geometric parameters (m; L) 
the ergodic components of the system have positive measure and on each of them 
the standard billiard ball flow has the K-property. 

Note 1. In general, if r ^ Ro, then Int M decomposes into a finite number 
of connected components. Our results extend to these cases, too. 

Note 2. As it will be seen in Section 2 (see Lemma 2.1) our system is 
isomorphic to a semi- dispersing billiard, and, what is more, is a semi-dispersing 
billiard in a weakly generalized sense, and, consequently, we can use the theory 
of semi-dispersing billiards. As has been proved in recent manuscripts by N. 
I. Chernov and C. Haskell, [C-H(1996)] on one hand, and by D. Ornstein and 
B. Weiss, [0-W(1998)] on the other hand, the K-mixing property of a semi- 
dispersing billiard flow on a positive ergodic component actually implies its 
Bernoulli property, as well. 

Note 3. As to the basic results (and history) concerning the Boltzmann 
and Sinai hypotheses we refer to the recent survey [Sz(1996)]. 

The basic notion in the theory of semi-dispersing billiards is that of the 
sufficiency of a phase point or, equivalently, of its orbit. The conceptual im- 
portance of sufficiency can be explained as follows (for a technical introduction 
and our prerequisites, see Section 2): In a suitably small neighbourhood of a 
(typical) phase point of a dispersing billiard the system is hyperbolic; i.e. its 
relevant Lyapunov exponents are not zero. For a semi-dispersing billiard the 
same property is guaranteed for sufficient points only! Physically speaking, a 
phase point is sufficient if its trajectory encounters in its history all possible 
degrees of freedom of the system. 

Our proof of the full hyperbolicity is the first major step in the basic 
strategy for establishing global ergodicity of semi-dispersing billiards as was 
suggested in our series of works with A. Kramli (see the references), and, in 
fact, some elements of that approach are also used in this paper. In the sense of 
our strategy for proving global ergodicity, initiated in [K-S-Sz(1991)] and ex- 
plained in the introductions of [K-S-Sz(1992)] and of [Sim(1992)-I] , there are 
two fundamental parts in the demonstration of hyperbolicity, once a combina- 
torial property, called richness of the symbolic collision sequence of a trajectory, 
had been suitably defined: 

(1) The "geometric- algebraic considerations" on the codimension of the man- 
ifolds describing the nonsufficient trajectory segments with a combinato- 
rially rich symbolic collision structure; 

(2) Proof of the fact that the set of phase points with a combinatorially non- 
rich collision sequence has measure zero. 
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In contrast to earlier proofs of ergodic properties of hard ball systems our 
method for obtaining full hyperbolicity will not be inductive; more precisely, we 
do not use the hyperbolicity of smaller systems. This has become possible since 
our Theorem 5.1 (settling step (2) here) which is a variant of the so-called weak 
ball- avoiding theorems, can now be proven without any inductive assumption. 
Nevertheless, the proof of our crucial Key Lemma 4.1, which copes with part 
(1), is inductive and, indeed, one of the main reasons to use varying masses is 
just that by choosing the mass of one particle to be equal to zero permits us 
to use an inductive assumption on the smaller system (but about its algebraic 
behaviour, only). As just mentioned, to prove our theorem we have introduced 
varying masses, and can only claim hyperbolicity for typical (N + l)-tuples 
(mi, ■ ■ ■ ,rriN',L) G IR^ +1 of outer geometric parameters. Furthermore, we will, 
in Section 3, complexify the dynamics. The complexification is required, on one 
hand, by the fact that our arguments in Section 4 use some algebraic tools that 
assume the ground field to be algebraically closed. Another advantage of the 
complexification is that, in the inductive derivation of Key Lemma 4.1, one 
does not have to worry about the sheer existence of orbit segments with a 
prescribed symbolic collision sequence: they do exist, thanks to the algebraic 
closedness of the complex field. 

The paper is organized as follows: Section 2 is devoted to prerequisites. 
Section 3 then describes the complexified dynamics, while in Section 4 we 
establish Key Lemma 4.1 — both for the complexified and the real dynamics 
- settling, in particular, part (1) of the strategy. Section 5 provides the 
demonstration of our Main Theorem through the aforementioned Theorem 
5.1, and, finally, Section 6 contains some comments and remarks. 

2. Prerequisites 
Semi-dispersing billiards and hard ball systems 

Our approach is based on a simple observation: 

Lemma 2.1. The standard billiard ball flow with mass vector m G is 
isomorphic to a semi- dispersing billiard. 

Because of the importance of the statement we sketch the proof. 

Proof. Introduce new coordinates in the phase space as follows: For (Q, V) 
G M let qi = ^Jrniqi, i>i = yjrnivi, where thus qi G T V L and J2iLi = 1- 
Moreover, the vector V = (v\,... ,vn) necessarily belongs to the hyperplane 
YliLi V™iVi = 0. Denote Y^L = YliLiiV™^ T l) and 

(2.2) Cij = \Q={ qi ,...,q N )e T^ L : 



fffU 



< 2r 
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and Q = T^x, \ Ui<j<j<jvC'jj. The equivalence relation over ^ is as 
follows: Q ~ > j, Q* if and only if there exists an a € T K such that for every 

i G [1, iV], = qi+yjmia allows us to define Q = Q/*m- Let Mrn,L = QxS d_1 
be the unit tangent bundle of Q, and denote by dp the probability measure 
const ■ dQ -dV , where <iV" is the surface measure on the (d— 1)— sphere § d_1 , and 
dQ is the Lebesgue- measure on frn,L- Now the standard billiard ball system 
(M, S*, n)m,L with mass vector m and the billiard system (M^l, <S*, /t) are 
isomorphic. 

Indeed, we can reduce the question to the case of one-dimensional parti- 
cles because for both models the velocity components perpendicular to the nor- 
mal of impact remain unchanged. The claimed isomorphy for one-dimensional 
particles, however, is well-known, and for its simple proof we can refer, for 
instance, to Section 4 of Chapter 5 in [C-F-S(1981)]. Thus, the point is that in 
the isomorphic flow (wirfi,,L, S*, fij , the velocity transformations at collisions 
become orthogonal reflections across the tangent hyperplane of the boundary 
of Q; see also (3.9) for the mentioned velocity transformation. 

The fact that the billiard system (M^l, <§ k , p) is semi-dispersing, is ob- 
vious, because the scattering bodies in Q are cylinders built on ellipsoid bases 
(therefore this system is a cylindric billiard as introduced in [Sz(1993)]). □ 

For convenience and brevity, we will throughout use the concepts and 
notation, related to semi-dispersing billiards and hard ball systems, of the 
papers [K-S-Sz(1990)] and [Sim(1992)-I-II], respectively, and will only point 
out where and how different masses play a role. 

Remark 2.3. By slightly generalizing the notion of semi-dispersing billiards 
with allowing Riemannian metrics different from the usual ones (see (2.2) of 
[K-S-Sz(1990)]), we could have immediately identified the standard billiard ball 
system (M, S*, with mass vector m with a (generalized) semi-dispersing 

billiard. (In the coordinates of the proof of Lemma 2.1, the Riemannian metric 
is 

N 

(dp) 2 = Y,(\\dK\\ 2 + \\ds\\ 2 ) 
i=i 

where ||ds|| 2 is the square of the natural Riemannian metric on the unit sphere 
S^ 1 .) Then the standard billiard ball system (M, S R , n)rh,L itself will be a 
(generalized) semi-dispersing billiard. The advantage is that — as is easy to 
see — the results of [Ch-S(1987)] and [K-S-Sz(1990)] remain valid for this class, 
too. For simplifying our exposition therefore, we will omit the change of coor- 
dinates of Lemma 2.1 and will be using the notion of semi-dispersing billiards 
in this slightly more general sense in which the model of the introduction is a 
semi-dispersing billiard. 
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An often used abbreviation is the shorthand S^'^x for the trajectory 
segment {S t x : a < t < b}. The natural projections from M onto its factor 
spaces are denoted, as usual, by it : M — > Q and p : M —> s N ' v ~ v ~ 1 or, 
sometimes, we simply write ir(x) = Q(x) = Q and p(x) = V(x) = V for 
x = (Q, V) G M. Any t £ [a, b] with S l x G dM is called a collision moment 
or collision time. 

As pointed out in previous works on billiards, the dynamics can only be 
defined for trajectories where the moments of collisions do not accumulate in 
any finite time interval (cf. Condition 2.1 of [K-S-Sz(1990)]). An important 
consequence of Theorem 5.3 of [V(1979)] is that — for semi-dispersing billiards 
under the nondegeneracy condition mentioned before our Main Theorem - 
there are no trajectories at all with a finite accumulation point of collision 
moments (see also [G(1981)] and [B-F-K(1998)]). 

As a result, for an arbitrary nonsingular orbit segment S^'^x of the stan- 
dard billiard ball flow, there is a uniquely defined maximal sequence a < t\ < 
ti < • • • < t n < b : n > of collision times and a uniquely defined se- 
quence o\ < o"2 < • • • < o n of "colliding pairs"; i.e. Ok = {ik,jk} whenever 
Q{t k ) = Tr(S tk x) e dC ikdk . The sequence S := Y,(S^x) := \a u a 2 , . . . ,a n ) is 
called the symbolic collision sequence of the trajectory segment S^'^x. 

Definition 2.4. We say that the symbolic collision sequence S = (o"i, . . . , a n ) 
is connected if the collision graph of this sequence: 

Or. ■= (V = {1,2, ... ,JV},f E := {{i k , j k } : where a k = {i k , j k }, 1 < k < n}) 
is connected. 

Definition 2.5. We say that the symbolic collision sequence £ = (<7i, ... ,a n ) 
is C-rich, C being a natural number, if it can be decomposed into at least C 
consecutive, disjoint collision sequences in such a way that each of them is 
connected. 

Neutral subspaces, advance and sufficiency 

Consider a nonsingular trajectory segment S^'^x. Suppose that a and b 
are not moments of collision. Before defining the neutral linear space of this 
trajectory segment, we note that the tangent space of the configuration space 
Q at interior points can be identified with the common linear space 

(2.6) Z = ^(w 1 ,w 2 ,...,w N ) e (R") N : ^m^ = oj . 
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Definition 2.7. The neutral space Mo(S^- a ' b ^x) of the trajectory segment 
S^'^x at time zero (a < < b) is defined by the following formula: 

Mv{S [a > h] x) = {W eZ: 3(5 > 0) such that Va G (-<5, <5), 

p (5 a (Q(x) + aW, V(x))) = p(S a x) and 

p(s b {Q{x) + uW,V{x))) =p(S b x)}. 

It is known (see (3) in Section 3 of [S-Ch (1987)]) that M (S^x) is a 
linear subspace of Z, and V(x) G Mo(S^ a ' b ^ x) . The neutral space Af t (S [a ' b] x) of 
the segment S^ a,b ^x at time t G [a, 6] is defined as follows: 

(2.8) Ai t {S^x) = AT (s^M . 

It is clear that the neutral space can be canonically identified with 

Afo(S^- a ' b ^x) by the usual identification of the tangent spaces of Q along the 
trajectory S^ -00 ' 00 ^ (see, for instance, Section 2 of [K-S-Sz(1990)] ). 

Our next definition is that of the advance. Consider a nonsingular orbit 
segment S^ a ' b ^x with symbolic collision sequence S = (a±, . . . , cr n ) (n > 1) as at 
the beginning of the present section. For x = (Q, V) £ M and W G Z, \\W\\ 
sufficiently small, denote T W (Q, V) := (Q + W, V). 

Definition 2.9. For any 1 < k < n and t G [a, b], the advance 

a(a k ) : Af t (S^x) R 

is the unique linear extension of the linear functional defined in a sufficiently 
small neighbourhood of the origin of Af t (S^ x) in the following way: 

a(a k )(W) := t k {x) - tkiS^Tw^x). 

It is now time to bring up the basic notion of sufficiency of a trajectory 
(segment). This is the most important necessary condition for the proof of 
the fundamental theorem for semi-dispersing billiards; see Condition (ii) of 
Theorem 3.6 and Definition 2.12 in [K-S-Sz(1990)] . 

Definition 2.10. 

(1) The nonsingular trajectory segment S^'^x (a and b are not supposed to 
be moments of collision) is said to be sufficient if and only if the dimension 
of Af t (S^x) (t G [a, b]) is minimal, i.e. dim AT t (S^x) = 1. 

(2) The trajectory segment S^'^x containing exactly one singularity is said 
to be sufficient if and only if both branches of this trajectory segment are 
sufficient. 

For the notion of trajectory branches see, for example, the end of Section 2 
in [Sim(1992)-I]. 
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Definition 2.11. The phase point x G M with at most one singularity is 
said to be sufficient if and only if its whole trajectory S^ -00 ' 00 ^ is sufficient, 
which means, by definition, that some of its bounded segments S^x are 
sufficient. 

In the case of an orbit S^ -00 ' 00 ^ with exactly one singularity, sufficiency 
requires that both branches of S^~°°'°°^x be sufficient. 

Connecting Path Formula for particles with different masses 

The Connecting Path Formula, abbreviated as CPF, was discovered for 
particles with identical masses in [Sim(1992)-II] . Its goal was to give an explicit 
description (by introducing a useful system of linear coordinates) of the neutral 
linear space Afo(S^ T '°^xo) in the language of the "advances" of the occurring 
collisions by using, as coefficients, linear expressions of the (pre-collision and 
post-collision) velocity differences of the colliding particles. Since it relied 
upon the conservation of the momentum, it has been natural to expect that 
the CPF can be generalized for particles with different masses as well. The case 
is, indeed, this, and next we give this generalization for particles with different 
masses. Since its structure is the same as that of the CPF for identical masses, 
our exposition follows closely the structure of [Sim(1992)-II] . 

Consider a phase point xq G M whose trajectory segment S^~ t '°^xq is 
not singular, T > 0. In the forthcoming discussion the phase point xq and the 
positive number T will be fixed. All the velocities, Vi{t) G W, i G {1, 2, . . . , N}, 
— T < £ < 0, appearing in the considerations are velocities of certain balls at 
specified moments t and always with the starting phase point xo («?(£) is the 
velocity of the i th ball at time t). We suppose that the moments and —T are 
not moments of collision. We label the balls by the natural numbers 1, 2, . . . , N 
(so the set {1, 2, ... , N} is always the vertex set of the collision graph) and we 
denote by ei, e2, . . . , e n the collisions of the trajectory segment S^~ t '°^xq (i.e. 
the edges of the collision graph) so that the time order of these collisions is just 
the opposite of the order given by the indices. More definitions and notation: 

1. ti = t(e{) denotes the time of the collision e,, so that > t\ > £2 > ■ ■ ■ > 

t n > ~T. 

2. If £ € R is not a moment of collision (— T < t < 0), then 

A ft (t) : M (S^x ) - R u 

is a linear mapping assigning to every element W G Mo(S^~ T '°^xo) the 
displacement of the i th ball at time £, provided that the configuration 
displacement at time zero is given by W. Originally, this linear mapping is 
only defined for vectors W G A/'o(S'[~ T ' ]xo) close enough to the origin, but 
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it can be uniquely extended to the whole space jVq(S^ T ' ^xq) by preserving 
linearity. 

3. a(ei) denotes the advance of the collision ef, thus 

a(e t ) : Af (S^ T ^x ) - R 
is a linear mapping {i = 1,2,..., n). 

4. The integers 1 = k(l) < k(2) < • • • < k(lo) < n are defined by the require- 
ment that for every j (1 < j < Iq) the graph {ei , e2 , . . . , e^j\ } consists 
of N — j connected components (on the vertex set {1,2,... , N}, as al- 
ways) while the graph {ei, e2, . . . , e^m-i} consists of N — j + 1 connected 
components and, moreover, we require that the number of connected com- 
ponents of the whole graph {e\, e2, . . . , e n } be equal to N — Iq. It is clear 
from this definition that the graph 

T = { e fc(l); e fc(2); • • • > e k(l )} 

does not contain any loop, especially lo < N — 1. 

Here we make two remarks commenting on the above notions. 

Remark 2.12. We often do not indicate the variable W 6 J\fo(S^~ T '°^xo) of 
the linear mappings Aqi(t) and a(ej), for we will not be dealing with specific 
neutral tangent vectors W but, instead, we think of W as a typical (run- 
ning) element of Afo(S^~ T '°^xo) and Aqi(t), a(e«) as linear mappings defined 
on J\fo(S^~ T '°^Xo) in order to obtain an appropriate description of the neutral 
space M (S^ T ^x ). 

Remark 2.13. If W € Af (S^ T ^x ) has the property A qi (0)[W] = Xvi{0) 
for some A £ R and for alH G {1,2,..., N} (here Uj(0) is the velocity of the i th 
ball at time zero), then a(efc)[W] = A for all k = 1,2, ... ,n. This particular 
W corresponds to the direction of the flow. In the sequel we shall often refer 
to this remark. 

Let us fix two distinct balls a, uj G {1,2,... , N} that are in the same 
connected component of the collision graph Q n = {ei,ez, . . . ,e n }. The CPF 
expresses the relative displacement Aq a (0) — Aq^O) in terms of the advances 
a(ej) and the relative velocities occurring at these collisions ej. In order to 
be able to formulate the CPF we need to define some graph-theoretic notions 
concerning the pair of vertices (a, uj). 

Definition 2.14. Since the graph T = {e^), e^), . . . , &k(i a )} contains no 
loop and the vertices a, uj belong to the same connected component of T, 
there is a unique path U(a,uj) = {/i, /2, . . . , fh} in the graph T connecting 
the vertices a and uj. The edges /j G T (i = 1, 2, . . . , h) are listed successively 
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along this path U(a,uj) starting from a and ending at uj. The vertices of the 
path H(a,u>) are denoted by a = B$,B\,B2, ...,Bf l = u indexed along this 
path going from a to u, so the edge fi connects the vertices -Bj_i and B{ 
(i = l,2,...,h). 

When trying to compute Aq a (0) — Ag a; (0) by using the advances a(ej) 
and the relative velocities at these collisions, it turns out that not only the 
collisions (i = 1, 2, . . . , h) make an impact on Aq a (0) — A(/ w (0), but some 
other adjacent edges too. This motivates the following definition: 

Definition 2.15. Let i € {1, 2, . . . , h — 1} be an integer. We define the set 
of Ai adjacent edges at the vertex Bi as follows: 

A = {ej : j € {1, 2, . . . , n} and (%) - ■ (tfo) - t(/ i+ i)) < and 

Bi is a vertex of ej } . 

We adopt a similar definition of the sets Ao, Ah of adjacent edges at the 
vertices Bq and Bh, respectively: 

Definition 2.16. 

.Ao = { e j '■ 1 < i < n an d i(ej) > and Bo is a vertex of e^}; 
^ = {ej : 1 < j < n and t(ej) > t(fh) and Bh is a vertex of e^}. 

We note that the sets Ao, Ai, . . . , Ah are not necessarily mutually disjoint. 

Finally, we need to define the "contribution" of the collision ej to Aq a (0) — 
Ag w (0) which is composed from the relative velocities just before and after the 
moment t(ej) of the collision ej. 

Definition 2.17. For i G {1,2, ... ,h} the contribution r(/j) of the edge 
/i G n(a,cj) is given by the formula 

v^Mf*)) ~ (*(/*)). if > and t(/ i+1 ) > 

1 



r(/i) = 4 



m Bl _ x + m Bi 



mB t - 1 + m B , 



[mB^ (v B ^ (t(fi)) ~ v Bi (t(/i))) + m Bi (t(/0) 
-«b 4 (*(/<)))]> if < < 

■K^t.W/')) -«£(*(/*))) (*(/*)) 

-«b 4 (*(/<)))]> if *(/<-i) < < *(/<+i). 
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Here (£(/?)) denotes the velocity of the Bf 1 particle just before the 
collision fi (occurring at time t(fi)) and, similarly, v~g.(t(fi)) is the velocity 
of the same particle just after the mentioned collision. We also note that, by 
convention, t(fo) = > and t(fh+i) = > t(fh). Apparently, the time 

order plays an important role in this definition. 

Definition 2.18. For i £ {0, 1,2, ... ,h} the contribution Ti(ej) of an edge 
ej £ Ai is defined as follows: 



Ti{ej) = S ign(t(fi)-t(f i+1 ))- 



m c 



m Bi + m c 

■ [K(%)) " v^(t( ej ))) - (v Bi (t(e 3 )) - vc(t(ej)))] 
where C is the vertex of ej different from i?,. 

Here again we adopt the convention t(/o) = > t(ej) (ej £ Ao) and 
t(fh+i) = > t(ej) (ej € Ah)- We note that, by the definition of the set 
Ai, exactly one of the two possibilities < t(ej) < t(fi) and t(fi) < 

t(ej) < occurs. The subscript i of T is only needed because an edge 

ej £ Ai 1 nAi 2 (i\ < 12) has two contributions at the vertices B^ and Bi 2 which 
are just the endpoints of ej. 

We are now in the position of formulating the Connecting Path Formula: 

Proposition 2.19. With the definitions and notation above, the following 
sum is an expression for Ag a (0) — Ag w (0) in terms of the advances and relative 
velocities of collisions: 

h h 

Ag Q (0)-A^(0)=^a(/ i )r(/ i ) + ^ «( e i) r i( e i)- 

i=l i=0 ej&Ai 

The proof of the proposition follows the proof of Simanyi's CPF (Lemma 
2.9 [Sim(1992)-II] ) with the only difference that Lemma 2.8 of [Sim(1992)-II] is 
replaced here by the following: 

Lemma 2.20. If e is a collision at time t between the particles B and C, 

then 

- = - -sw) + - >SW) 

and 

t,+ (t) - v~(t) = m C [(t,+(t) - «+(*)) - (v B (t) -vc(t))]. 
iHB ~r 'fie 

Proof. The lemma is an easy consequence of the conservation of momen- 
tum for the collision e. □ 
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Remark 2.21. In Section 3, we will complexify the system, and also allow 
complex masses, in particular. It is easy to see that the CPF of Proposition 
2.19 and the whole discussion of this subsection will still hold for complex 
masses if we assume, in addition, that for every e = {B, C}, occurring in the 
symbolic collision sequence of S^-~ T '°^xo, mgl mc / 0. 



3. The complexified billiard map 
Partial linearity of the dynamics 



The aim of this section is to understand properly the algebraic relationship 
between the kinetic data of the billiard flow measured at different moments of 
time. 

From now on we shall investigate orbit segments S^°' T ^xq (T > 0) of the 
standard billiard ball flow ^M, {S l },fi, m, . We note here that M = Q x £, 

where the configuration space Q is as defined right after (1.1) and £ is the 
velocity sphere 

N N } 



£ = { (Vl, . . .,V N ) <= 



Y^rriiVi = and ^^mj||fj|| 2 = 1 > 
i=i i=i J 

introduced in Section 1. Also note that in the geometric-algebraic considera- 
tions of the upcoming sections we do not use the equivalence relation of the 
introduction. This is why we work in M rather than in M. Later on even the 
conditions YliLi m i v i = and YliLi m ill v ill 2 = 1 wm be- dropped (cf. Remark 
3.14). 

The symbolic collision sequence of S[ ' T lxo is denoted by 



T] 



£o) = (0l,---,0"n)- 



The symbol vf = vf{tk) G W denotes the velocity qi of the i th ball right 
after the k th collision = {ik,jk} (1 < ik < 3k < N) occurring at time 
k = 1, . . . , n. Of course, there is no need to deal with the velocities right before 
collisions, since v~(t k ) = v?(t k -i) = u* -1 . As usual, q\ G T U L = W/L ■ If 
denotes the position of the center of the i th ball at the moment t k of the k th 
collision. 

Let us now fix the symbolic collision sequence S = (a±, . . . , a n ) (n > 1), 
and explore the algebraic relationship between the data 



and 



{ q t\vt 1 \i = l,...,N} 
{gf,^|z = l,...,iv}, 
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k = 1, ... ,n. By definition, the data {qf , vf\i = 1, . . . , iV} correspond to the 
initial (noncollision) phase point xq. We also set to = 0. 

Since we would like to carry out arithmetic operations on these data, 
the periodic positions qf G are not suitable for this purpose. Therefore, 
instead of studying the genuine orbit segments S^°' T ^xq, we will deal with their 
Euclidean liftings. 

Proposition 3.1. Let S^'^xq = {%(*), Vi(t)\ <t <T} be an orbit seg- 
ment as above, and assume that certain pre-images (Euclidean liftings) <ji(0) = 
q® G W of the positions %(0) G T U L = W /L ■ Z u are given. Then there is a 
uniquely defined, continuous, Euclidean lifting {q~i(t) G W \ < t < T} of the 
given orbit segment that is an extension of the initial lifting i = 1, . . . , iV}. 
Moreover, for every collision o k there exists a uniquely defined integer vector 
a k eZ" - named the adjustment vector of o k - such that 



(3.2) 



t ik -q k jk -L-a k 



4r 2 = 0. 



The orbit segment Co = [qf, vf \ i = 1, . . . , N; k = 0, . . . , n} is called a lifted 
orbit segment with the system of adjustment vectors A = (a±, . . . , a n ) G 7L nv . 

The proof of this proposition is straightforward and we omit it. 

The next result establishes the already-mentioned polynomial relation- 
ships between the kinetic data 

{^- 1 ,^- 1 |i = i > ...,jv} 

and 

{qtvf\i = l,...,N}, 

k = 1, . . . , n . 

Proposition 3.3. Using all notation and notions from above, one has 
the following polynomial relations between the kinetic data at o k ~\ and o~ k , 
k = 1, . . . ,n. In order to simplify the notation, these equations are written as 
if o k were {1,2}: 



(3-4) v? = vr 1 , t*{l,2}, 

(3.5) m\v i + m 2 v 2 = miv^ 1 + m 2 v\~ x 

(conservation of the momentum), 
(3.6) 

„,k „fc _ „fc-l „fc-l ! /„fe-l ,,k-l~k ~k T n \ ( ' ~k ~k j n \ 
V\~v 2 —v x -v 2 - -v 2 , q x - q 2 - L ■ a k J ■ \q x - q 2 - L ■ a k J 
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(reflection of the relative velocity determined by the elastic collision), 
(3-7) tf = $- 1 + Tkv!- 1 , i = l,...,N, 

where the time slot = — t^-i (to := 0) is determined by the quadratic 
equation 

(3.8) l^- 1 - g*" 1 + r k (t,*" 1 ~ v*- 1 ) ~L-a k f = 4r 2 . 

The proof of this proposition is also obvious. 

Remark. Observe that the new velocities v\ and v\ can be computed from 
(3.5)-(3.6) as follows: 

(oq\ „fc _ „,fc-l m 2 /„,fc-l „,fc-l.~fc ~k j n 

(3-9) v l =v 1 - ■ r < v 1 -v 2 ■q 1 - q 2 - L- a k 

v 2~v 2 + 2r2(mi + m2) { v i "-2 ,?i-fe-^-a fc 
Qi — Q2 — L ■ ak) ■ 



We also note that the equations above even extend analytically to the case 
when one of the two masses, say mi, is equal to zero: 

(3.10) 



- k ~ 1 1 /^k- 1 „ k ~ 1 -n k n k T n \ ( n k n k T n \ 

Vi -v 2 , q l - q 2 - L ■ at ) ■ I q± - q 2 ~ L> • afc I , 



2r 2 



We call the attention of the reader to the fact that in our understanding 
the symbol ( . ; . ) denotes the Euclidean inner product of ^-dimensional real 
vectors and || . || is the corresponding norm. 

The complexification of the billiard map 

Given the pair (£, A) = (ci, . . . , a n ; a±, . . . , a n ), the equations (3.4)-(3.8) 
make it possible to iteratively compute the kinetic data {qf, \ i = 1, . . . , iV} 

by using the preceding data jgf -1 , v*~ l \i = l,...,iV"j. Throughout these 
computations one only uses the field operations and square roots. (The latter 
one is merely used when computing as the root of the quadratic equation 
(3.8). Obviously, since - at the moment - we are dealing with genuine, real 
orbit segments, in this dynamically realistic situation the equations (3.8) have 
two distinct, positive real roots, and is the smaller one.) Therefore, the data 
{gf , vf | i = 1, . . . , A/"} can be expressed as certain algebraic functions of the 
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initial data vf\ i = 1, . . . , iV"}, and the arising algebraic functions merely 
contain field operations, square roots, the radius r of the balls, the size L > 
of the torus and, finally, the masses m; as constants. 

Since these algebraic functions make full sense over the complex field C 
and, after all, our proof of the theorem requires the complexification, we are 
now going to complexify the whole system by considering the kinetic variables, 
the size L, and the masses as complex ones and by retaining the polynomial 
equations (3.4)— (3.8). However, due to the ambiguity of selecting a root of 
(3.8) out of the two, it proves to be important to explore first the algebraic 
frame of the relations (3.4)-(3.8) connecting the studied variables. This is what 
we do now. 

The field extension associated with the pair (£,«4) 

To avoid misunderstanding we immediately stress that the field exten- 
sions K = K(T,;A) to be defined below will also depend on a sequence r = 
(to, . . . , t„_i) of field elements to be introduced successively in Definition 3.11. 
If we also want to emphasize the dependence of K on r, then we will write 
K = K(Z;A) = (Z,A,t). 

We are going to define the commutative function field K = K (£; A) gen- 
erated by all functions 

{(Qi)j, (vf)j\i = l,...,N; k = 0,...,n; j = 

of the lifted orbit segments corresponding to the given parameters 

(£,*4) = (oi, . . . ,(J n ;ai, . . . ,a n ) 

in such a way that the field IK (£; A) incorporates all algebraic relations among 
these variables that are consequences of equations (3.4)-(3.8). (Here the sub- 
script j denotes the j th component of a ^-vector.) In our setup the ground 
field of allowed constants (coefficients) is, by definition, the complex field C. 
The precise definition of K n = K (£; A) is: 

Definition 3.11. For n = the field Kq = K(0;0) is the transcendental 
extension C(B) of the coefficient field C by the algebraically independent formal 
variables 

B = (v%, m u L\i = l,...,N;j = l,...,iy}. 

Suppose now that n > and the commutative field K n _i = K (£'; A') has 
already been defined, where £' = (a±, . . . , cr n _i), A' = (a±, . . . , a n -±). Then 
consider the quadratic equation b n r^ + c n T n + d n = in (3.8) with k = n as a 
polynomial equation defining a new field element r n to be adjoined to the field 
K n _i = K(T,';A'). (Recall that the coefficients b n , c n , d n come from the field 

Kn-l.) 
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There are two possibilities: 

(i) The quadratic polynomial b n x 2 + c n x + d n is reducible over the field K n _i. 
Then K n = K ra _i. 

(ii) The polynomial b n x 2 + c n x + d n is irreducible over the field K n _i. Then 
we define K n = K (£; A) as the extension of K n _i = K (£'; A') by the root 
T n of this irreducible polynomial. 

Remark 3.12. The importance of the field K n is underscored by the fact 
that this field encodes all algebraic relations among the kinetic data that follow 
from the polynomial equations (3.4)-(3.8). Furthermore, Proposition 3.3 gives 
an iterative computation rule for successively obtaining the kinetic variables 



k = 0, . . . ,n. The field K n = K(E,„4) is, after all, the algebraic frame of 
such computations. We note that - everywhere in what follows - the sym- 
bols ( . ; . ) and || . || 2 do not refer to an Hermitian inner product but, rather, 
(x; y) = Ylj=i x jVj an d IM| 2 = Sj=i x p so that these expressions retain their 
polynomial form. 

Remark 3.13. Later on it will be necessary to express each kinetic variable 
(<2^) . and (w^) . as an algebraic function of the initial variables of B. This raises, 
however, an important question: Which one of the two roots T/% of (3.8) should 
be considered during these computations? This is no problem if the polynomial 
bkX 2 + CkX + dk, defining r^, is irreducible over the field because the two 

roots of this polynomial are algebraically equivalent over K^-i, and any of them 
can be chosen as r^. However, if the polynomial bkX 2 + CkX + dk is reducible 
over the field IfC^-i, then it is necessary to make a decision and assign to one 
of the two roots (both in the field now) the role of r^. This is what we 

do. The result is the field lK n = K(E;^4) endowed with a distinguished n-tuple 
(ri, . . . ,T n ) of its elements. The field K n with the distinguished n-tuple r is 
denoted by K ra = !€(£; A; f). Then the algebraic object K(S; A; f) completely 
controls the whole process of the iterative computation of the kinetic variables. 

Note that in the genuine, real case the two roots of bkX 2 + CkX + dk are 
distinct, positive real numbers, and - by the nature of the billiard dynamics - 
the chosen value is always the smaller one. 

Finally, we mention here a simple case when the polynomial bkX 2 +CkX+dk 
is reducible over the field Namely, this takes place whenever Ok-\ = Cfc 

and cifc-i = cifc. In that case the value Tk = is clearly a solution of the above 
polynomial bkX 2 + c^x + dk-, and, being so, that polynomial is reducible over 
the field see also 3.31-3.33. 




Vi\l = 



l,...,Ar}, 
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Remark 3.14. The reader may wonder why we did not postulate the al- 
gebraic dependencies YliLi m i v i = 0> YliLi m i\\ v i l| 2 = 1- The answer is the 
following: The definition of the neutral linear space together with its charac- 
terization via the Connecting Path Formula (CPF, see Section 2 of this paper 
or Lemma 2.9 of [Sim(1992)-II]), i.e. the partial linearity of the billiard flow, 
are invariant under 

(1) all uniform velocity translations (adding the same vector to all velocities) 
and 

(2) time rescalings. 

Remark 3.15. All fields occurring in this paper are only defined up to an 
isomorphism over the coefficient field C. Therefore, the statements like "the 
field K2 is an extension of Ki" should be understood as follows: K\ is a subfield 
of K2 after the natural identification of the generating variables bearing the 
same name. For the necessary notions, properties, and results from the theory 
of field extensions and Galois theory, the reader may look up, for instance, the 
book by I. Stewart, [St(1973)]. 

Remark 3.16. As said before (see Definition 3.11), the collection B of 
elements of the field K(S; A; f) is a base of transcendence in that field over the 
subfield C, and the degree of the extension K(S; A; f) : C(£>) = K n : Kq is a 
power of two. Moreover, each of the following sets is a generator for the field 

I,-- - ,N; j = 1,.. . ,i/}U{tl, . . . ,r n }; 
l,...,N; j = l,...,u; k = 0,...,n}. 

Remark 3.17. The above procedure of constructing field extensions is 
closely related to the classical theory of geometric constructions by a ruler and 
compass, and this is not surprising: the billiard trajectory is constructed by 
intersecting a straight line with a sphere and then reflecting it across the tan- 
gent hyperplane of the sphere. This is a sort of classical geometric construction 
in v dimensions. For the details see, for instance, Sections 57 and 60 of [VDW 
(1955)], or Chapter 5 of [St(1973)]. 

Let us fix the pair (S; A) and the n-tuple r = (t±, . . . , r n ) of elements of 
the field K(S,^4). We are now defining a (2u + l)N + 1-dimensional complex 
analytic manifold & = Q,(T>,A, f) and certain holomorphic functions qf, vf 
: n -> C u , mi : Cl -> C, L : Cl -> C, T k : Q -> C, i = 1, . . . , TV, k = 0, . . . , n 
(tq is not defined). But first we introduce: 



(a) CUBU{Ti,...,T n }; 

(b) CU{mi, . . .,m N ,L}U{(^)j, {v?)j\i = 

(c) CU{mi, . . . ,m N , L}um )j, (v^\i = 
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Definition 3.18. Define the domain D = D {Z,A,t) C C^ 2u+1 ^ n+1 as the 
set of all complex (2v + 1)N + 1-tuples 

( v i)j, mi, L\i = l,...,N; j = l,...,u} 

for which 

2 



<*-** 



(a) The leading coefficient 

ratic equation (3.8) is never zero, k 
equivalent to 



and the discriminant of the quad- 
1, . . . , n. (The latter condition is 



1 - v k . 1 ; q k 



L-a k )^0 



Ik 3k 

see also Remark 3.29.) And 

(b) rrii k + nij k ^ for k = 1, . . . , n, provided that the iterative computation 
of the kinetic variables is carried out, 

K fc )il * = 1, ■ ■ ■ ,N; j = 1, . . • , v; k = 0, . . . ,n} . 

as governed by the polynomial equations (3.4)-(3.8) and by the fixed se- 
lection of elements f. 

Remark 1. When inverting the dynamics, we see that 



and 



tk 



3 k 



4 



■jr 1 



4^ 



L ■ a k 



Remark 2. Note that we require the validity of (a) for any branch of 
the square root function, when computing the r^'s with irreducible defining 
polynomials (3.8). 

Lemma 3.19. The complement setC^ 2u+1 ^ N+1 \D is a proper, closed, al- 
gebraic subset of C,( 2l/ + 1 ) N + 1 j especially a finite union of complex analytic sub- 
manifolds with codimension at least one, so that the open set D C C( 2u+ V N+1 
is connected and dense. 

Proof. We have the polynomial equations (3.4)-(3.8) making it possible to 
set up an algorithm for iteratively computing the kinetic variables measured 
at different times. The point is not just this iterative algorithm, but also 
its invertibility (time reversibility). The inverse process has similar algebraic 
properties, for it just means time reversal. One can easily prove by an induction 
on the length n of S that the iteratively defined dynamics determines a several- 
to-several mapping with maximum rank between a nonempty Zariski open set 
of the kinetic variables with superscript 1 and a nonempty Zariski open set of 
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the kinetic variables measured right after the last reflection (with superscript 
n). (A Zariski open set is the complement of a closed algebraic set defined 
as the simultaneous zero set of finitely many polynomials.) Then one can go 
ahead with the induction from n to n + 1, because all possible obstructions 
to extending the process only occur on proper algebraic submanifolds of the 
Zariski open set of kinetic variables with superscript n. □ 



Definition 3.20. We define Q, (T,,A,t) as the set of all complex (2v(n + 
l)N + N + n + l)-tuples 



for which these coordinates 

(1) are interrelated by the equations (3.4)-(3.8), 

(2) respect the choices of Tk prescribed by f (cf. 3.11, and also 3.13), 

and the vector 



of initial data belongs to the set D (T,,A, f) defined above. 

It is clear now that qf, v\ : & — > C v , rrii : Q, — > C, L : Q, — > C, 
and Tk : £1 — ► C (i = 1, . . . , N, k = 0, . . . , n) are holomorphic functions on 
0, = Cl(E,A, f). The complex analytic manifold 0, endowed with the above 
holomorphic functions can justifiably be considered as the complexification of 
the (S, A, r)-iterated billiard map. 

Remark. The careful reader has certainly noted the following principal 
observation: the r^'s figure in two different roles. They first denote succes- 
sively chosen field elements of the extensions introduced in Definition 3.11, 
and secondly they also denote complex- valued functions on (l (or multivalued 
complex functions on D; cf. Definitions 3.20 and 3.18). 

Remark. It is very likely that the manifold £1 is connected. Nevertheless, 
in our proof we do not need this connectedness, and are not going to pursue 
the goal of proving it. 



Fix a base point uj G £1 (S, A, r). The tangent space 7^0 of Q at u consists 
of all complex vectors 



u = 




...,n;j = l,...,u;l = l,...,n 




The complex neutral space jV(uj) 



x = (<5i 



q°,SvlSm i ,5L\i = l,...,N)eC^ N + 1 ; 
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thus can be naturally identified with the complex vector space 
Set 
(3.21) 

= Lc= (5$, 5v°, Smi, 6L\i = 1,. . . ,N) €%Q 

D x (vf) = D x (mi) = D X {L) = for i = 1, . . . , TV; k = 0, . . . , n 

where D x (.) denotes the directional differentiation in the direction of x. 
Clearly, the complex linear subspace M(uj) of c( 2u + 1 ) N + 1 [ s the proper complex 
analogue of the neutral space of a genuine, real orbit segment. 

Remark. Since 8v® , 8rrii , and SL must be equal to zero for a neutral vector 
x G M(co), we shall simply write x = (Sq^, • • • , o~q%) instead of indicating the 
zero entries. 

The proof of the following proposition is completely analogous to the one 
for the real case, and therefore we omit it. 

PROPOSITION 3.22. For every tangent vector x G J\f(u) (u e & (T,,A,f)) 
and 1 < k < n, 

(3.23) D x (g*- 1 - = (a k (x) - a k ^(x)) • (v^ 1 - v^ 1 ) , 

(3.24) D x (t k ) = -a k (x), 
and 

(3.25) D x (r k ) = D x (t k - t k -i) = a k -i(x) - a k (x), 

where, by definition, to = = ao(x). The functions a k : M(uj) — > C (k = 
1, . . . , n) are linear Junctionals. The name of the functional a k is the advance 
of the k th collision a k = {i k ,j k }; see also Section 2. 

This result markedly shows that, indeed, the vector space M(uj) describes 
the linearity of the (S, A, r)-iterated billiard map. 

Remark 3.26. It is a matter of simple computation to convince ourselves 
that all assertions in Section 2 of the present paper and in [Sim(1992)-II] per- 
taining to the neutral space Af(u>) remain valid for the complexified dynamics. 
Here we only point out the two most important statements from those results: 

(1) The vector space M(to) measures the ambiguity in determining the or- 
bit segment to purely by its velocity history {v£ (u) : i = 1, . . . , N; k = 0, . . . , n} 
and the outer geometric data (m, L). In other words, this means that if, locally 
in the phase space Q. (in a small open set), two phase points u>\ and 002 have 
the same velocities {v^{uj-i) = vf(u>2), i = 1,...,N; k = 0, . . . , n) and geomet- 
ric parameters (m, L), then the initial data of u)\ and ll> 2 can only differ by a 
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spatial translation by a vector from the space M(uji) = N(uJ2)- This statement 
is obviously reversible. 

(2) The reflection laws for the neutral vectors are exactly the same as for 
velocities. 

The Connecting Path Formula (CPF, Lemma 2.9 of [Sim(1992)-II], or 
Proposition 2.19 in this paper) is applicable, and it enables us to compute the 
neutral space N(uo) via solving a homogeneous system of linear equations (now 
over the field C) 

n 

(3.27) ^2a k T ik = 0, i = 1, . . . ,n + Ps - N, 

k=l 

where each equation in (3.27) is a z^-dimensional complex vector equation, the 
coefficients = Tik(uj) G C u are certain linear combinations of relative ve- 
locities vf(ui) — Vj(u>) (the coefficients of those linear combinations are just 
fractional linear expressions of the masses, see Definitions 2.17-2.18) and, fi- 
nally, P-£ denotes the number of connected components of the collision graph 
of S. For a more detailed explanation of (3.27), see the beginning part of 
Remark 4.6 in the next section. 

Remark 3.28. We note here that the ordering of the moments of colli- 
sions t(fi) (which plays a significant role in Definitions 2.17-2.18) is no longer 
meaningful over the unordered field C. However, the use of the corresponding 
inequalities in 2.17-2.18 is purely technical/notational. Those inequalities only 
serve to introduce the combinatorial ordering of the collisions fi (which are cr's) 
given by the indexes of the symbolic collision sequence £ = (o"i, 02, . . . , cr n ). 

Remark 3.29. As a straightforward consequence of equation (3.6) in Propo- 
sition 3.3, we obtain that v i ~ 1 = v\ for i = 1, . . . , N (i.e. the k th reflection 
does not change the compound velocity) if and only if 

«- 1 -«i- 1 ;i-*- L - a *> = °. 

and this is just the case of a tangential collision. An easy calculation shows 
that the mentioned tangentiality occurs if and only if the two roots of (3.8) 
coincide. 

It is obvious that dimcAA(w) is at least v + 1 (as long as not all velocities 
are the same), because the flow direction and the uniform spatial translations 
are necessarily contained by J\f(uj). 

Definition 3.30 (Sufficiency). The orbit segment to £ (l(T,,A, f) is said 
to be sufficient if and only if dimcAA(w) = v + 1. 

Finally, the main result of the next section (Key Lemma 4.1) will use the 
following: 
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Definition 3.31. The triple (E,A,t) has Property (A) if and only if the 
following assertion holds: 

For every pair of indices 1 < k < I < n for which = o~i, = ai, and 
(jj PI cife = for j = k + 1, . . . , I — 1, we have that 77 ^ — ^i=fc+i r i i n the 
algebraic selection r of the time slots. (Note that in this case the quadratic 
polynomial b\x 2 + c\x + d\ is reducible over the field because the time slot 
r\ = — Yl l f=k+i T 3 G ^-1 i s automatically a solution of (3.8).) 

Remark 3.32. It is obvious that the schemes (Y,,A,t) corresponding to 
genuine, real billiard trajectory segments enjoy Property (A). Indeed, in those 
cases even the equations a/~ = o~i, = ai, and o~j PI a/~ = (J = k + 1, . . . , I — 1) 
cannot hold simultaneously. Concerning the necessity for introducing Property 
(A), see Example I at the end of the next section. 

Remark 3.33. For combinatorial schemes satisfying Property (A), the 
quadratic polynomial bix 2 + cix + di is reducible over the field (i.e. = 
Ki) only in the case described above in 3.31. This will be a consequence of 
Main Lemma 4.21 (cf. Corollary 4.36). 

4. The notion of richness. Richness implies sufficiency 

The result of this section is 

Key Lemma 4.1. There exists a positive number C(N) (depending merely 
on the number of balls N > 2) with the following property: If a symbolic se- 
quence S = (<7i, . . . , a n ) is C(N)-rich and the triple (E, A, f) has Property (A) 
(see Definition 3.31), then dimcAA(w) = u + 1 (i.e. the trajectory is sufficient) 
for almost every to £ Q(T,,A,t). The real version of this result is also valid. 
(For the definition of a C-rich symbolic sequence, see 2.5.) 

Remark 4.1/a. Plainly, the exceptional zero measure subset of 17 (T,,A,f) 
mentioned in this lemma must actually be a countable union of proper, analytic 
submanifolds. 

The aim of Sections 3 and 4, in particular that of Key Lemma 4.1, is to 
settle part (1) of the strategy detailed. The methods are mainly algebraic, and 
for an easier separation of our dynamical vs. algebraic arguments, we are now 
going to describe what are the only properties of the dynamics (and nothing 
more) that we use. 

(1) The dynamical properties of an orbit segment S^'^xo with a symbolic 
collision sequence E are completely characterized by equations (3.4)-(3.8) if we, 
in addition, specify t/% as the smaller of the two, necessarily positive, roots of 
equation (3.8). In the inductive steps of Lemmas 4.9, 4.39, and 4.40 we also use 
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the remark to Proposition 3.3 saying that the collision laws (3.9), themselves 
consequences of (3.5)-(3.6), extend analytically to the case of zero mass; cf. 
(3.10). 

(2) Since our goal is to show that orbits with a rich collision structure 
are typically sufficient or, equivalently, their neutral subspaces have typically 
minimal possible dimension, we also need an algebraic characterization of the 
neutral subspace. This is attained through the CPF (see Proposition 2.19 and 
the system of linear equations (3.27)), and our arguments also use some simple 
inferences about the CPF and the neutral subspaces from [Sim(1992)-II] . 

(3) As we have seen in Section 3, the algebraic formalisms of points (1) 
and (2) can be complexified, and our forthcoming discussion is based upon this 
complexification. 

Remark 4.1/b. In the course of proving the key lemma we will see that the 
thresholds C(N) (N > 2) can actually be defined by the recursion C(2) = 1, 

C(N) = f • max{C(iV - 1),3} (N > 3); that is, C(2) = 1, C(N) = -^-j 

for N > 3. We note that, as proved in [K-S-Sz(1991)], for A = 3 the value 
C(3) = 2 already does the job. 

Proof of Key Lemma 4.1. All the rest of this section is devoted to the 
inductive proof of the key lemma, and which be split into several lemmas. 

First of all, we prove the original, complex version of the key lemma by 
an induction on the number TV. Once this is completed, we will show that the 
validity of the key lemma carries over to the real case. 



The dichotomy. Algebraic characterization of sufficiency 

Lemma 4.2. Consider the open, connected and dense subset D of 
C( 2u + 1 ) N + 1 formed by all possible complex (2u + 1)N + 1-tuples 

x(lo) = ((g?(w)) . , (v^(oo)) j , rm(u), L{u)\ i = 1, . . . , N; j = 1, . . . , v) (wefi) 

of initial data of the complex (S, A, t) -dynamics; see Definitions 3.18 and 3.20. 
We claim that there are finitely many complex polynomials P\ (x) , . . . , P s (x) 
(x £ ^ 2u+1 ^ N+1 ), canonically determined by the discrete parameters (Y,,A,t), 
such that for every x £ D (£, A, f) the following implication holds true: 

(4.3) 

(3 to £ U (S, A, f) such that x = x(uj) and dimcAA(w) > v + 1^ 

=> Pi(x) = P 2 (x) = ■■■ = P s (x) = 0. 



HARD BALL SYSTEMS ARE COMPLETELY HYPERBOLIC 



59 



Remark. The left-hand side in the above equivalence precisely says that 
some branch of the multiple-valued (E, A, r)-dynamics with initial data x is 
not sufficient, see also Remark 2 to Definition 3.18. 

Proof. This proof is a typical application of the Connecting Path Formula 
(see Proposition 2.19 here or Lemma 2.9 and Proposition 3.4 of [Sim(1992)-II]). 
The latter asserts that 

(4.4) dim c AA(u;) = v ■ P s + dim c {ai, . . . , a n }, 

where, as said before, Ps denotes the number of connected components of the 
collision graph of E, and {a±, . . . , a n } is shorthand for the complex linear space 
of all possible n-tuples (ai(x), . . . , a n (x)) of advances, x £ M(uj). 

Remark 4.5. It is worth noting here that our present formula (4.4) dif- 
fers from its counterpart in Proposition 3.4 of [Sim(1992)-II] by an additional 
term v. This is, however, due to the fact that in the present approach we 
no longer have the reduction equation ~Y^a=\ m i^Qi = 0- It follows easily from 
(4.4) that the sufficiency is equivalent to dimc{o!i, • • • , ct n } = 1; i.e., suffi- 
ciency means that all advances are equal to the same functional. Note that 
dime {or, • • • , a n } = 1 obviously implies Pe = 1. 

A simple, but important consequence of the Connecting Path Formula is 
that the linear space {cti, . . . , a n } is the solution set of a homogeneous system 
of linear equations (3.27). 

Remark 4.6. Note that, when applying the mentioned lemma, the left- 
hand side of the Connecting Path Formula has to be written as the relative 
velocity of the colliding balls multiplied by the advance of that collision. It fol- 
lows immediately from the exposition of [Sim(1992)-II] that equations of this 
sort (arising from all CPF's) are the only constraints on the advance function- 
al The reason why this is true is that the fulfillment of all CPF's precisely 
means that the relative displacement (variation of position) of every pair of 
particles right before their new collision is parallel to the relative incoming 
velocity of these particles (see also (3.23)), and this fact guarantees that the 
variation of the newly formed relative outgoing velocity will also be zero, just 
as required in (3.21). 

Also note that the number n + Ps — N of equations in (3.27) is equal to 
the number of all collisions = {ik,jk} for which and jk are in the same 
connected component of the collision graph of {<7i, 02, . . . , 

The determinant D(M) of every minor M of the coefficient matrix of 
(3.27) is a homogeneous velocity polynomial and, therefore, D(M) is a holo- 
morphic function on £1 (E, A, f). Thus, there are finitely many velocity polyno- 
mials (actually, determinants) Ri(to),... ,R s (uj) such that their simultaneous 
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vanishing is equivalent to the nonsufficiency of the orbit segment u. By us- 
ing the equations (3.4)-(3.8) (which define the complex dynamics recursively), 
each of these velocity polynomials R%(uj) can be - essentially uniquely - written 
as an algebraic function /j (x(ui)) = fi(x) of the initial data x, and these alge- 
braic functions fi(x) only contain (finitely many) field operations and square 
roots. By using the canonical method of successive elimination of the square 
roots from the equation fi(x) = 0, one obtains a complex polynomial Pi(x) 
such that 

Pi{x) = <^> fi(x) = 0. 

We emphasize here that 

(1) we understand the relation fi{x) = in such a way that fi{x) becomes 
zero on some branch of the square root function when evaluating the 
multiple-valued algebraic function fi(x); cf. Remark 2 to Definition 3.18; 

(2) the equivalence is claimed to only hold on D; cf. the formulation of Lemma 
4.2. 

The process of eliminating the square roots from an equation f(x) = 
(of the above type), however, requires a bit of clarification. The algebraic 
function / has a natural representing element / in the field K n (see Definition 
3.11). Suppose that / G K k but / ^ K k _ ± . Then the field element / 6 K k 
can be written uniquely in the form / = a + 6\/<5, where a, b, 5 G b ^ 0, 

\f5 and 

M^'-^C'-C 1 -^ 2 



< x -< 1 



2 



2 



-4r 2 



is the discriminant of the quadratic equation (3.8) with the unknown r k . 

In the first step of the successive elimination of the square roots from 
the equation f(x) = 0, we switch from the field element f = a + by/6 6E IfCfc 
(representing /) to the new field element 

/' = (a + bV6^j ■ (a - bVS^j = a 2 - b 2 5 G K*_i, 

which is a representing element of an algebraic function f'(x) containing square 
roots than f(x). Corresponding to the switching / h> /', we also replace the 
original equation f(x) = by the new equation f'(x) = 0. Then we continue 
this process until we arrive at a field element 

m = ^gi = c(f), 
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which is a complex rational function — of the initial variables x. Of course, 

here we assume that the polynomials P(x) and Q(x) have no nontrivial com- 
mon divisor. 

Actually, as directly follows from Proposition 3.3 and from the method 
of eliminating the square roots (presented above), in each denominator Qi(x) 
(1 < I < s) associated to the algebraic function fi(x) there are only factors of 
the type m; + rrij or \\vf — Vj\\ 2 , l<i<j<N,0<k<n. Therefore, none of 
the denominators Qi{x) vanishes on D. This explains property (2) formulated 
above. 

As it also follows from the method of eliminating the square roots and from 
the actual form of the equations (3.4-3.8), the coefficients of the polynomials 
Pi(x) are real. 

By using some simple facts from Galois theory, we can easily detect an inti- 
mate relationship between the rational function g(x) £ Ko just constructed and 
the product a of all conjugate elements of / £ K n . (The conjugates of / (over 
the base field Ko) are sitting in the normal hull ~K n : K.q of the field extension 
K n : Kq. The element a 6 Kq is obviously the constant term in the normalized 
minimal polynomial of / over the field Kq. For the elementary concepts and 
facts from Galois theory, the reader is referred to the book [St(1973)].) By 
establishing a natural matching between all conjugates of a 2 — b 2 5 and the 
conjugates of the ordered pair (a + bVS, a — bs/d), we can easily see that there 
exists an integer I such that g = a 2 . (In case I < this should be understood 
as a = (g) 2 .) Especially, a and g have the same irreducible factors both in 
the numerator and in the denominator. Another consequence is that the field 
element / is zero if and only if P{x) = 0. 

So far we have seen that for each velocity determinant Ri(u>) (i = 1, . . . , s) 
one can find a canonically determined polynomial Pi{x) with the following 
property: For every vector x £ D (£, A, f) the equation Pi(x) = is equivalent 
to the equation fi(x) = on some branch of the square root function when 
evaluating the algebraic function fi(x). Hence Lemma 4.2 follows. □ 



Corollary 4.7 (dichotomy corollary). As another consequence of the 
above identity g = a 2 ' , we obtain that if P = 0, then f{x) = for every x € 

C (2zH-l)AT+l and 

on every branch of the square root function when evaluating the 
multiple-valued algebraic function f(x). This also means that if all polynomials 
Pi, . . . ,P S in Lemma 4.2 happen to be identically zero, then no orbit segment 
to G (l(T,,A,f) is sufficient. 

Therefore, thanks to the algebraic feature of the dynamics, there is a di- 
chotomy: Either every (S, A, r)- orbit segment is nonsufficient, or almost every 
such trajectory segment is sufficient. 
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Remark. We note that the reason why we do not have an equivalence in 
(4.3) is that - for a typical x G c( 2jy+1 ) 7V+1 - the branches, on which fi(x) = 0, 
may be different for different values of i. 

According to the above lemma, in order to prove Key Lemma 4.1, it is 
enough to show that at least one of the polynomials Pi(x) is nonzero. This is 
exactly what we are going to do by an induction on the number N > 2. 

Lemma 4.8. Pertaining to the polynomials Pi(x) associated with Ri(uj) = 
fi (x(u))) in Lemma 4.2, the following holds: 

If the algebraic function fi(x) is identically zero on a nonempty, open sub- 
set B ofMP v+1 ^ N+1 , then fi is algebraically trivial, i.e., fi = or, equivalently, 
Pi(x) = 0. 

Proof. Denote the minimal polynomial of fi over the base field Ko = C(x) 
by M{x) = J2j=o c j x ^ c j G ^0- Suppose, on the contrary, that fi vanishes on 
an open ball B of R^+^+i. Since f. fuimis the ident ity ^=0^(0 = 0, 
the vanishing of fi on the ball B implies that the element cq of the rational 
function field M,q vanishes on B, as well. Therefore, Co = in IKo, so that M(x) 
is divisible by x. The irreducibility of M yields M(x) = x, i.e., fi = 0. Hence 
the lemma follows. □ 

The substitution {m^ = 0}. Derived schemes 

The next lemma describes the (£, A, f)-orbit segments in the case when 
the iV th particle is "infinitely light" compared to the others, i.e., when mjy = 
but mi ■ ... ■ mN-i / 0. 

Lemma 4.9. Suppose that = and m\ • . . . • mjy_i / 0. Assume 
further that there exist indices 1 < p < q < n such that N £ a p , N G a q and 
N ^ Oj for j = p + 1, . . . , q — 1. Then an orbit segment u £ Q (X, A, r) is 
sufficient (i.e. dimcAA(u;) = v + 1) if 

(1) the {1,2,... ,N -\}-part 

{qf(Lo), vf(u), mi(w), L(w) | i = 1, . . . , N - 1; k G /jv} 

of the orbit segment is sufficient as an orbit segment of the particles 
1, . . . , N — 1 and 

(2) the relative velocities v p n (uj) — v i ^(uj) andv q N ~ 1 (uj)—v c !~ 1 (Lu) are not parallel. 

(Here, as usual, i p (i q ) is the index of the ball colliding with the iV th 
particle at a p (cr q ), and the index set In C {0,1, ... ,n} contains and those 
indices i > for which iV ^ crj.) 
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Remark. It is worth noting that, as follows easily from the conditions, 
v p N (iv) = v q ^ l {uj). 

Proof. First of all, we observe that the motion of the -/V th ball with zero 
mass has absolutely no effect on the evolution of the {1, . . . , N — l}-part of the 
trajectory segment, and this statement is also valid for the time evolution of the 
{1, . . . , iV-l}-part {5q^\ i = 1, . . . , N - 1} of a neutral vector (Sq^, ... , 5q%) G 
M(lo); see also the equations (3.4), (3.7)-(3.8) and (3.10). 

Consider now an arbitrary neutral vector 5Q = (Sql, . . . , 5q%) £ Af(uj). 
According to (1) and the above principle, we can modify the neutral vector SQ 
by a scalar multiple of the flow direction and by a uniform spatial displacement 
in such a way that Sq® = for i = 1, . . . , N — 1, and then Sqf = remains true 
for the whole orbit segment, i = 1,. . . ,N — 1, k = 0, . . . , n. The neutrality of 
SQ with respect to the entire orbit segment means, however, that 

(4.10) Sq p N = a p ■ (v P N (u) - uf p (w)) and Sq"" 1 = a q ■ (u^V) ~ "'"V)) ■ 

According to our hypothesis (2), the consequence of (4.10) and of the obvious 
equation Sq p N = Sq^ 1 is that a v = a q = 0, and thus 5Q = 0. Hence the lemma 
follows. □ 

The crucial part of the inductive proof of Key Lemma 4.1, i.e. the sub- 
stitution m^v = (and L = 0, too), requires some preparatory thoughts and 
lemmas. The formulations given for the substitution tun = can be easily 
adapted to the substitution L = 0. Assume that a combinatorial-algebraic 
scheme (£, A, f) is given for the iV-ball system {1,2,..., iV} so that Property 
(A) holds (see Definition 3.31). In the proof of Key Lemma 4.1 we want to use 
some combinatorial-algebraic (N — l)-schemes (T,',A',t') that 

(i) govern the time evolution of the {1, 2, ... , iV— l}-part of some (T,,A, f)- 
orbit segments with an infinitely light -/V th ball, i.e. mjv = and 

(ii) enjoy Property (A). 

In the construction of such (N — l)-schemes (£' ,A' ,t'), the so-called de- 
rived schemes, we want to follow the guiding principles below, which also serve 
as the definition of the derived schemes: 

Definition 4.11. We say that the (N — l)-scheme (X',„4,',r') is a scheme 
derived from (T,,A,t) by putting mjy = if it is obtained as follows: 

(1) First of all, we discard all symbols cij, dj, and Tj from (T,,A,t) that 
correspond to a aj containing the label N. This also means that we retain 
the other symbols dj, aj without change: they only get re-indexed, due to the 
dropping of the other symbols; 

(2) As far as the selection of r' is concerned, we want to retain all algebraic 
relations among the variables that are encoded in the original scheme (£, A, r), 
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i.e. all algebraic relations that follow from the scheme (T,,A,t) and from 
tun = 0; 

(3) We want the derived scheme (T i ',A',t / ) to enjoy Property (A); see 



It is straightforward that, following just the instructions in (l)-(2) above, 
one can easily construct such schemes (T,,A,t) which, perhaps, do not have 
Property (A). 

Remark 4.12. Part (2) of the above definition means the following: Ini- 
tially, the algebra of (T,,A,r) consists 

(1) of the variables 



(2) of the equations (3.4-3.8), and finally 

(3) of the choices of the signs of the square roots in the solutions of the 
quadratic equations (3.8) whenever the equations are reducible (cf. Re- 
mark 3.13). 

To obtain the algebra of the derived scheme (£', A', r'), we put first 
tun = and then t[ = Y^j=k+i T i whenever for some 1 < k < I < n we 
have N £ <Tk U <ti, N € Dj'=fc+i ° j (here the indexing of the time slots t[ of 
(E',A',t') is lacunary). After trivial substitutions, the quadratic equation for 
t i leads to the natural quadratic equation for Tj , and then, finally, we can can- 
cel all the variables {(q%)j, (v%)j, m^l < k < n, 1 < j < v}, and all the 
equations they figure in as well. The collection of the remaining variables and 
equations with the corresponding vector r' make (£', A', ?). 

Remark 4.13. Below we show the inheritance of the reducibility = 

¥Li =^ IFCJ = Kj. We note, however, that this result is not used in the proof of 

Key Lemma 4.1. We merely include it - along with its nice consequence: part 
(3) of Lemma 4.37 - in order to give the reader a broader view of the subject. 
The subsequent proof of Lemma 4.40 only uses the weaker part (2) of 4.37 
which, in turn, does not need the above mentioned inheritance of reducibility. 

Another interpretation of part (2) of the above definition is the following 
one: Once we have a complex rational expression g(z\, . . . ,z p ) of arbitrary 
elements 



so that the same expression of these variables is zero in the field K n = 
K(T,, A,t), then it must also be zero in the field K' n . An important case 



3.31. 




z 1 ,...,z p £K' n = K{Z',A',T > ) 



HARD BALL SYSTEMS ARE COMPLETELY HYPERBOLIC 



65 



(when (2) has special significance) is when N £ crjt U 07 (1 < k < I < n), 
N e nj=fe+i a ji an d EQ-i = Kj. Then the time slot t{ = Y?j=k+i T i (here the 
indexing of the time slots t! of (T,',A',f') is lacunary.) should be expressed 
as a rational function of the earlier defined (£, A, r)-variables (the variables 
constituting the field UQ_i), just as it is indicated by the rational dependence 
of 77 on the variables constituting the field Kj_i. Therefore, r[ = Y?j=k+i T i ls 
a rational expression 

(4.14) t[ = g (x;q° N ,v° N ,m N ;T 1 , . . . ,ti^) 

of the initial variables (x; q%, v%, ttin) € c( 2v + 1 ) N + 1 an d the time slots 
Ti, . . . , ti-\. Note that the identity (4.14) is also true if we impose additionally 
the side condition = 0: 

(4.15) 7-; = g (x ; q° N , v° N ,0; rf, . . . , r^) . 

However, the {1, . . . , N— l}-dynamics evolves independently of the data of the 
N th ball with zero mass tun- This means that 

»(f;gSr(l) > <(l) > 0;7f(l) > ...,7f_ 1 (l)) 

= < 7 (f;4(2),^(2),0;r 1 (2),...,r z °_ 1 (2)) 

whenever all cumulative sums Y^=k!+i r j(l) are ecma l to the correspond- 
ing cumulative sums YljLkx+i r j(2) for 1 < fci < A;2 < Z, iV ^ 07^ U a^, 
N G fljtfa+i ®j • Thi s precisely means that the rational function gyx;q^,v%,0; 

T l > • • • ) T Z-lJ 

(i) does not contain the variables (q%) ^, (v%) ^ (1 < /3 < u), and 

(ii) only depends on the above cumulative sums YljLki+i rj 1 = r^j. 

(The indexing of t'® is still lacunary!) Thus, we get that g(x; q° N , v%, 0; r°, . . . , 
r°_ 1 ) is, in fact, a rational expression of x and the time slots preceding r^' . 
This shows that r^ G i.e. = follows from the similar reducibility 

Kj-i = Kj. 

The main point is the following one: Once 77 is the root of a reducible 
polynomial, it is selected "artificially", following the instructions encoded in 
the scheme f. Now we want to avoid the selection of the other root when 
defining the derived scheme (£',,A / ,r / ). This is why (2) is a requirement, and 
not a corollary, and this is why, in Remark 4.12, we use the wording "the time 
slot t[ should be expressed." 

The fact that the reducibility of the polynomial defining 77 is inherited dur- 
ing a substitution tun = followed from a general principle: If some algebraic 
function in Kj_i (the discriminant of the quadratic polynomial defining 77) is 
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the square of another such algebraic expression, then this square relation re- 
mains even "more true" if one annihilates a variable, i.e. considers that square 
relation on a smaller set. It is now the nature of our model that mjy = im- 
mediately kills all variables involving the motion of the N th particle, since the 
considered new discriminant of the quadratic equation defining r[ is obviously- 
independent of those variables. 

Remark 4.16. (a) Note that the construction of such schemes (E' , A' ,7'), 
by just requiring (1)— (2), is far from unique: It is very possible that, originally, 
for some k we have \Kk : Kfc-i| = 2 (i.e. the polynomial (3.8) determining 
Tfc is irreducible), but after the substitution = the corresponding field 
extension collapses to a trivial, degree-one extension. (The paramount example 
for this phenomenon is just the set-up of the upcoming Main Lemma 4.21.) 
If this phenomenon of "reducible splitting" takes place, then we have to make 
our choice: Which one of the two (now distinguishable) roots of the reducible 
(3.8) to include in t"! 

The requirement (A) for the derived scheme (£', A', ?'), however, ties our 
hands when selecting this root in the scenario of Main Lemma 4.21, and it 
actually implies the uniqueness of the derived scheme (£', A', t 1 ) (see Corollary 
4.36). 

The phenomenon of this reducible splitting can also be well represented by 
the following, simple paradigm: Let a € C be a complex parameter (playing the 
role of rriTv). Consider the Riemann-leaf of the function y = f(z) = \J z 2 + a, 
i.e. the algebraic variety defined by the equation y 2 — z 2 — a = 0. This is 
irreducible (undecomposable) for a / 0, and it splits into two components 
(y — z = and y + z = 0) for a = 0. However, given a sequence (ik — > 
(a& ^ 0), every point (z, y) on the curve y 2 — z 2 = turns out to be a limit of 
a sequence (z^, yu) with y\ — z 2 — = 0. 

(b) The paradigm just given also makes it possible to illustrate why (2) in 
Definition 4.11 is a requirement and not a corollary. Indeed, consider the field 
Ko = C(z,mjv); i.e., ~Kq is the complex field extended by two transcendental 
elements z and m^. Let Ki := Ko(ri) = Ko( \J z 2 + 1 + mjy) be the exten- 
sion of Ko by the solution t\ of the irreducible equation t 2 — z 2 — 1 — mjy 
= 0. Let us now extend Ki by the root T2 of the quadratic polynomial 
x 2 — 5t\x + 6t 2 . This polynomial is reducible over the field K4: its two roots 
are 1t\ and 3t\. Now the extended field K2 will be just equal to Ki, and 
the selection of a designated root (as T2) of the above polynomial amounts to 
selecting either 2t\ or 3t\ . Observe that the root of the considered reducible 
polynomial could have been an arbitrary rational expression of the elements 
z, and T\. In a general reducible step, such a rational expression shows the 
unique way of computing the actual r„ from the initial data and the Tj's with 
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i = 1, . . . ,n — 1. However, when assigning actual complex values to the field 
element t\ = \/ z 2 + 1 + mjv, we have two possiblities. Summing everything 
up, our selection does not mean giving values to these r's but, rather, telling 
if we take (formally) 2t\ or 3t\ as ti- (This selection determines the discrete 
algebraic structure, which becomes transparent if we think about repeated cr's 
with identical adjustment vectors.) When defining the derived scheme this 
means that once we have chosen 3i/ z 2 + 1 + tun (instead of 2\/ z 2 + 1 + tun) 
we still want to choose 3V z 2 + 1 after the substitution mjy = 0, and do not 
want to switch suddenly to the other root 2a/ z 2 + 1. 

Trans versalities of the degeneracies 

Definition-notation 4.17. Denote by <& p (1 < p < n) the degeneracy 

2 ' / 

ff p — v?^ =0 and by ^ p (1 < p < n) the degeneracy defined by (vf p — v p ^ ; 

Qi p ~ Q P j v ~ L " a pj = 0- (These degeneracies are represented here with post- 
collision velocities, but in view of Remark 1 to Definition 3.18 they also have 
an analogous representation in terms of pre-collision velocities.) 

Remark 4.18. (a) In the following arguments we will use a slightly different 
working domain D = D p+ i )Tt _i, which is not exactly the one of Definition 
3.18. Instead, in general, the working domain D p+ i ira _i describes the orbit 
segments that are nondegenerate from the (p+ l)-st collision up to the (n — 1)- 
st one (inclusive). Here p < n — 1, and if p = n — 1, then, by definition, 
.Dp_|_i in _i = c( 2!/ + 1 )^ v + 1 , xhe reference time (the upper index of the initial 
data) is now p + 1; i.e. the outgoing kinetic variables are taken at a p+ \. 

(b) Observe that if for some p < p' 

(4.19) o p = oy and, for p < j < p', Uj PI a p = either a p or 0, 

then by Remark 1 to Definition 3.18 $ p = <Ey. Similarly, if for some p < p' 
(4.20) 

(dp, a p ) = (dp/, a p i) and forp < j < p either (oj,aj) = (a p , a p ) or Oj n o p = 0, 

then necessarily ^ p = VEy. 

(c) In the upcoming Main Lemma transversality statements for <£ p (or 
typ) and <!>„ (or fy n ) will be claimed on D p+ i >n -i. Because of our previous 
observations, these assertions can be nonvacuous only if the symbolic pair 
(St, A) is minimal in the sense that 

(1) if <£p is involved in the statement, then there is no p' : p < p' < n such 
that (4.19) holds (this requirement is called -minimality from the left); 

(2) if fy p is involved in the statement, then there is no p' : p < p' < n such 
that (4.20) holds (this requirement is called ^-minimality from the left); 
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Analogous requirements can be formulated when <E> n (* n ) is involved, and will 
be called <3> (or -minimality from the right. 

In order to show that there really exist derived schemes (that is, fulfilling 
not only (l)-(2), but also (3) of Definition 4.11), we need to prove the following: 

Main Lemma 4.21. Assume that (T,,A, f) enjoys Property (A). We claim 

that 

(I) The codimension-one submanifolds defined either by rrii = (1 < 
i < N) or by L = can not locally coincide with any degeneracy <E> P or ^ p , 
1 < p < n; 

(II) On D p+1 . n — i, the degeneracy ^> p cannot locally coincide with the de- 
generacy ^> n (1 < p < n — 1) whenever ((<t p , . . . , a n ), (a p , . . . , a n )) is ^-minimal 
from the left and $ -minimal from the right; 

(III) On D p+ i jra _i, the degeneracy ^ p cannot locally coincide with the 
degeneracy whenever {{a p , . . . , a n ), {a p , . . . , a n )) is ^-minimal from both 
sides; 

(IV) The degeneracies rrii = (1 < i < N), L = intersect D p+ i tn _i 
(1 < p < n — 1); moreover the degeneracies <3> n , ^> p and ^ n intersect 
D p+ i^i (1 < p < n — 1) whenever ((<r p , . . . , cr n ), (a p , . . . , a n )) is ^-minimal 
from the left or from the right, or is ^ -minimal from the left or from the right, 
respectively; 

(V) If (E,A, t) is ^ -minimal from the right, then there exists a (nonde- 
generate) orbit segment 

u = {<??H, vf(w), mi{u), L(u)\i = l,...,N; a = 0, 1, . . . , n - 1} 

G fi (<7i, . . . , a n -\; a\, . . . , a n -\; t\, . . . , r n _i) = f2 n _i 

for which the discriminant 



» - <-» 2 • ( ^V) - C» ~ L ■ a n 



4r 2 



o/ (3.8) (u>i£/j k = n) is equal to zero and rrii(uj) ^ (i = 1, . . . , iV), L(lo) ^ 0; 
(VI) If 1 < k < n, a k = o n , a k = a n , N a n , and 

(i) for k < j < n either aj n a n = or N £ Oj , 

(ii) £/iere exists a jo iratt k < jo < n and aj Q n <r n / 0, 

i/ien |K ra : K n _i| = 2 or, equivalently, the polynomial (3.8), defining r n , is 
irreducible over the field K n _i; 
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(VII) Ifl<k<n,(Jk = o n , a-k 7^ a n , and for every k < j < n, ajf]a n = 0, 
then again |K n : IK n _i| = 2. 

(VIII) On D p+ i jra _i, i/ie degeneracy Q p cannot locally coincide with the 
degeneracy $ n (1 < p < n — 1) whenever ((a p , . . . ,a n ), (a p , . . . ,a n )) is <J>- 
minimal from both sides. 

Remark 4.22. (a) Assertion (VI) of Main Lemma 4.21 is tailor-made to 
ensure that whenever the substitution tun = is carried out at every "critical" 
step with the possibility of selecting the wrong T' n (and thus violating Property 
(A)), we are not forced to do so, because at an irreducible r n either branch of 
the corresponding square root of the discriminant of (3.8) may be selected for 
continuing the dynamics. Quite similarly, assertion (VII) guarantees the same 
thing for the substitution L = 0. 

(b) The careful reader has certainly observed that statement (VIII) is 
separated from the quite similar assertions (II)— (III) of the Main Lemma. The 
reason is the following: By a slight improvement of the inductive arguments 
verifying the transitions N — 1 — ► N when establishing (II)— (III) and (VIII), 
one can strengthen these claims by only requiring in each of them minimality 
from the right. The advantage of these stronger statements is that then (VIII) 
is not needed in the inductive proofs of (II)-(IV). As a result, (V)-(VII) can 
be proved completely without (VIII). 

Proof. 1. First we prove (I) and for the sake of brevity, only for the 
case nii = 0. The handling of the submanifold L = is completely analogous. 
Assume that there is a counterexample, and p (the time of degeneracy) has the 
smallest possible value. Push forward the manifold rrii = up to the collision 
&p-i by the algebraic dynamics: the form of the equation rrii = defining this 
manifold does not change. (This push-forward can be made without difficulty 
by the minimality of the index p.) In this way we can reduce the proof of (I) 
to the case p = 1. However, the statement is obviously true in that case. □ 

For later use we prove the following: 

Corollary 4.23. If for the combinatorial scheme (E,„4, f) and for some 
ujq £ ^n-i, S(ujq) = 0, then there is an u>* £ (l n -i such that 5(u>*) = 0, none 
of the masses mi(u/), . . . ,toat(c<j*) is zero and, moreover, L(uj*) / 0. 

Proof. Indeed, the conditions imply that for a co dimension-one algebraic 
subset Z of D ^ n _i^ it is true that for every x £ Z there exists an w 6 & n -i 

such that x(u) = x and 6(u>) = 0. By (I), Z is not even locally identical with 
any submanifold {m 8 = 0} (i = 1, . . . , N). Consequently, we can select a small 
perturbation x* £ Z of xq = x(loo) and an u* £ £l n -i such that x(uj*) = x* , 
5(uj*) = 0, and none of the masses mi(w'), . . . ,rriN(u)*), nor L(u*), vanishes. □ 
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2. The implication (II) and (III) and (IV) => (V) is obvious since, by the 
argument used to prove (I) , it is no problem to avoid the additional degeneracies 
nii + rrij = 0, 1 < i < j < N. 

3. We prove next that the existence of an w G £l n -i with = 
(guaranteed in (V)) implies \K n : K n _i| = 2. This will immediately prove the 
implications (V) =^ (VI) and (V) =^ (VII). (It is easy to check that the 
conditions in (VI) or (VII) actually imply the precondition of (V).) 

Let, therefore, to G £l n -i be an orbit segment with 5{to) = 0. By re- 
labeling the balls, we can assume that a n = {1,2}. In a small neighborhood 
Uq C f^n-i °f w the discriminant <5(o/) from (V) (u/ G Uq) is an analytic 
function of uj' or, equivalently, S(u') is an analytic, algebraic function f(x) of 
the initial variables x = x{uo) G C<> +1)Ar+1 . 

Consider now the "normal vector of collision" 

w = w(u>) = q\{oj) — q% (<^>) — L • a n 

= - S~TV) -L-a n + r n M • [vrH") - vT 1 ^)] , 

where T n (oj) is just the double root of (3.8) with k = n. Execute now the small 
perturbations lu(e) G Uq, (e G C, |e| is small) as follows: 

ST 1 ( W (£)) = gT 1 H + e-w; 

q"' 1 (u(s)) = ft' 1 (u) , t = 2 JV; 

v?- 1 ("(e)) = kT 1 H , « = !,...,*. 

(Note: If 1 G o" n _i, then the configuration perturbation defined above ought to 
be modified by an appropriate scalar multiple of the flow direction in order to 

= 4r 2 .) Easy geometric consideration yields 



ensure 
that 



Zn-l ri n ~ 1 T n 



£=0 



That is, the holomorphic function g(e) := f (x (u)(e))) of £ G C (|e| is small) 
has the properties g(0) = 0, </(0) 7^ 0. 

Suppose now, indirectly, that the algebraic function f(x) G K n _i is a 
square (h(x)) 2 in K n _i. Consider the function z(e) := /i (x (lo(e))) of the single 
complex variable e G C, |e| small. It is clear that z(e) is a meromorphic function 
of e, being a rational function of the variables {qf,vf,rrii,L\i = 1, . . . , N; 
a < n}, and thus being the ratio of two holomorphic functions. This is the 
point where we use the indirect assumption h(x) G K n _i. However, the square 
[z(e)] 2 = g{e) of the meromorphic function z(e), obviously, cannot have the 
properties g(0) = 0, g'(0) / 0. 

This proves that, indeed, \K n : K n _i| = 2. □ 
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So far we have proved that (II), (III) and (IV) together imply assertions 
(V)-(VII) of Main Lemma 4.21. 

4. We prove now (II)-(IV) ((VIII) will be settled later in point 5.) Let us 
first observe that (II)-(IV) (and similarly (VIII)) are consequences of somewhat 
weaker statements 

(II*) On D2, n -i, the degeneracy <J>i cannot locally coincide with the de- 
generacy \l/ n whenever ((ci, . . . , cr n ), (a±, . . . , a n )) is ^-minimal from the left 
and ^-minimal from the right; 

(IIP) On L>2,n-i, the degeneracy \&i cannot locally coincide with the 
degeneracy ^ n whenever ((<ti, . . . , a n ), (a±, . . . , a n )) is ^-minimal from both 
sides; 

(IV*) The degeneracies rrii = (1 < i < N), L = intersect D2,n-ij 
moreover, the degeneracies $i, $ n , and \I/ n intersect L>2,n-i whenever 
((<7i, . . . , a n ), (ai, . . . , a n )) is ^-minimal from the left or from the right, or 
is ^-minimal from the left or from the right, respectively; 

(VIII*) On D2,n-i) the degeneracy <l?i cannot locally coincide with the 
degeneracy Q n whenever ((01, . . . , cr n ), (cti, . . . , a n )) is ^-minimal from both 
sides. 

Indeed, we only have to shift the index (time) of the initial variables from 
2 to p + 1 and use the observation: If the triple (S, A, r) has Property(A) and 
(T,',A') = (a p ,...,a n ; a p ,...,a n ) (1 < p < n) is an end segment of (E,A), 
then the S'-restriction of every (S, A, r)-orbit segment a; is automatically a 
(S', ^4', r')-orbit segment with some r', so that (S', .4', r') enjoys Property (A). 
Observe, moreover, that if rj (p < j < n) is reducible in the truncated dynamics 
(Hi 1 , A', t'), then Tj is certainly reducible in the original dynamics (T,,A, f). In 
other words, any possible additional irreducibility in the truncated dynamics 
permits a larger variety of orbit segments to study than just the restrictions 
of the orbit segments of (T,,A, f) to the collisions a p ,...,a n . Therefore it is 
sufficient to establish the claimed transversalities for (£', A', r'), and this is 
exactly the content of (II*)-(III*) and (VIII*). 

We will establish now (IP)-(IV*) (and (VIII*) later in point 5) by using 
a double induction: we assume that the statements are true for any number 
N' < N of balls and arbitrary length n' < n such that the pair (n',N') is 
not identical with (n,N). More precisely, when establishing (IV*) for the pair 
(n, N), we will only use (IV*) for pairs (n', N) with n' < n. On the other hand, 
when demonstrating (II*)— (III*) (and (VIII*)), the validity of (IF)-(IV*) (and 
also of (VIII*)) for pairs (n',N f ) with n' < n and N' < N will be exploited. 
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The proof is formulated in such a way that we take first the outer induction 
on n, and then the statements for this fixed n are shown by using an inner 
induction on N . 

Observe that, for an arbitrarily fixed number N of balls, assertion (IV*) is 
certainly true for the value n once (II*)-(IV*) and (VIII*) are known to hold 
for smaller values of n. Indeed, for the degeneracies rrtj = (1 < i < N) and 
L = 0, the statement is obvious. Our inductive hypothesis (i.e. the validity of 
(II*)— (III*) and (VIII*) for all smaller values of n) guarantees that the claimed 
intersections of the singularities with the working domain Z?2,n-i are nonempty, 
as long as the appropriate minimality conditions of (IV*) hold. 

Let us turn to the inductive proofs of (II*)— (III*) (and to that of (VIII*) 
later in point 5). 

The base of the induction on n is n = 2: This case is a straightforward 
inspection. 

Assume now that n > 3, and that (IP)-(IV*) and (VIII*) have been 
proved for all smaller values of n. The base of the induction on N is N = 
2 or 3 or 4, depending on the form of (£,„4), when establishing (II*) and (III*). 
In both proofs the inductive steps will first be presented, and afterwards the 
base cases will be treated. 

By re-labeling the balls, we can obtain that a n = {1,2}, o\ = {i,j}, 
1 < i < j < 4, and if j = 4, then i = 3 ("densely packed" indices). 

4/a: Proof of (II*) by using the inductive assumption, i.e. by assuming that 
(IP)-(IV*) and (VIII* ) are valid for pairs (n',N') satisfying n' <n,N' < N. 
We argue by contradiction. Suppose, therefore, that, on Z)2,n-i; $i locally 
coincides with the degeneracy \I> n defined by the equation S(lo) = 0, i.e. their 
intersection contains a nonempty, smooth, codimension-one submanifold of 
Di,n-\- We apply our favorite: 

Globalization/ substitution argument. In the sense of Remark 4.18, our 
variables (€ £>2,n-i) will now be represented at reference time 2, and the 
(algebraic) dynamics exists (strictly) between the collisions 1 and n. By using 
two quadratic equations: 

(1) the one expressing that o\ is a collision of two particular balls; 

(2) the other defining our considered o\ -degeneracy <I>i, 

we can eliminate two variables, say {qf) 1 and {vf) v by expressing them as 
some algebraic functions f\ (y) , fa (y) of the other variables 



a = l,...,N;l<(3<u^\{(qf) 1 , (^)J 
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measured at the collision o-i. The kinetic variables of y measured at 02 (he. 
with superscript 2) are now serving as the initial, or reference variables in the 
upcoming globalization/substitution argument. As usual, the algebraic func- 
tions fi(y) and f2(y) only contain field operations and square roots. Then, by 
using the recursive feature of the (S, A, f)-dynamics and the inductive assump- 
tion of 4.21, we can express the discriminant 5(u) as another (also essentially 
unique) algebraic function g(y) of the variables y. 

We argue by contradiction. Assume, therefore, that the degeneracies <I>i 
and fy n locally coincide in the working domain Di, n -\- The local coincidence 
of these degeneracies precisely means: 

(*) The algebraic expression g{y) is identically zero on some nonempty, open 
set of vectors y G C^ 1 ^' 1 . 

Then, according to Lemma 4.8, this provides, on L>2,n-i, the global contain- 
ment: $1 C ^ n . This statement remains valid after the substitution = 0, 
provided that it was valid in the original set-up. We remind the reader that 
the almost everywhere (nondegenerate) definability of the algebraic dynamics 
subjected to the side condition my = easily follows from (I). 

Remark 4.24. Since the substitutions tun = 0, L = play a basic role in 
the forthcoming arguments, it is very important to note here that, thanks to 
the inductive hypothesis of 4.21, we can assume that Property (A) is inherited 
by the algebraic billiard trajectory segment subjected to the side conditions 

= and L = 0. Indeed, as we mentioned in Remark 4.22, parts (VI)-(VII) 
are designed to ensure the inheritance of Property (A). 

Let us return now to our proof of the inductive step of (II*). Assume, of 
course, the inductive hypothesis. We have, moreover, already assumed indi- 
rectly that 3>i and ty n coincide locally in D2,n-i- By the inductive assumption, 
$1 and ty n intersect D2, n -i, whereas by the indirect assumption and by the 
globalization argument 

(4.25) $1 C t> n 

holds on D<i, n -\- If> under the substitution 777, jy = (say), the degeneracy <&! 
(and thus *$> n , too) intersects D' , the working domain (cf. Remark 4.20) for 
the derived scheme (T i ',A',t') of (E,A,t) corresponds to my = 0, then the 
contradiction to the inductive hypothesis (II*) is evident. For every k with 
a', G denote by & k and the projections of and ^ to the phase space 
of the (N — l)-ball system corresponding to taking m/y = and cancelling 
the data of the N th ball (as many times before, the indexing in the derived 
schemes is again lacunary). 

Consider now the case when, under the substitution = 0, the degen- 
eracy or does not intersect D' . The relation n D' = can only 
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hold if the symbolic sequence £' is not ^-minimal from the left and, similarly, 
^' n nD f = can only happen if the symbolic sequence (£', A') is not ^-minimal 
from the right. 

In any case, there exists a subsegment (£" , A" ,t") of (E', v 4 / ,r / ), which 
is minimal from both sides. Moreover, by denoting by p" and n" the indices 
of its initial and final elements (in the original indexing of (E,A,t)), then, 
on one hand, a[ = a' p „, and for every a'j G X/ with 1 < j < p", o'- n crj = 
either a[ or 0, and on the other hand, (a' n ,a' n ) = (a' n „,a' n „), and for every 
a'j £ S' with n" < j < n either (crj, a^) = (cr^, aJJ or n^ = 8. Note that, 
by virtue of Remark 4.18/b, the just formulated conditions also imply $' x = & p „ 
and %>' n „ = ty' n . Also, by inductive assumption (IV*), & pl , n ZV+i, n »-i / 
and ^' n „ n D p ii +l n ii_i ^ and further, by (II*), they cannot coincide locally 
in D p fr +l n ii_i. 

When completing the argument we will be exploiting the following lemma, 
which is a limiting variant of the observation of point (b) of Remark 4.18 and 
can be proved easily. 

Lemma 4.26. Consider a symbolic collision sequence (T,,A,f) which sat- 
isfies Property (A) and is ^-minimal from the left and ^ -minimal from the 
right. Suppose we set = and denote by (£',„4',t') the corresponding 
derived scheme. Now, there are two statements (whose duals for $ -minimality 
from the right and ^-minimality from the left are also true): 

(1) Assume that a[ = a' p , (= {1,2}, say) 1 < p' < n, and for every a'j £ £' 
with 1 < j < p' one has a'j n a[ = a[ or (lacunary indices). Take a 
sequence ujk £ & (£; A; r) such that 

(i) all kinetic data of oj k measured right after o\ have a limit, 

(ii) lim m N (uj k ) = 0, 

fc^oo 

(hi) lim rrniuk) / 0, 1 < i < N - 1, 

fc— *oo 

(iv) lim nii(uj k ) + mj(u k ) / 0, 1 < i < j < N - 1, 

fc— >oo 

(v) lim ||vi(wfc) - vl((jj k )\\ 2 = 0. 

fc— >oo 

Now limfc^oo (uj k ) - v\ (w fc )|| 2 = 0. 

(2) Assume that (o-' n ,,a' n ,) = (o-' n ,a' n ) (a' n = {1,2}, say) 1 < n' < n, and for 
every cr'- G £' with n' < j < n one either has that (o~j,a'j) = (o-' n ,a' n ) or 
a j H cr^ = is true. 

Take a sequence uj k G (l (T,;A;t) such that (i)-(iv) hold (with o\ replaced 
by a n in (i)), and further 

(v)' lim (v?(u k ) - v%(u k ); q?(cu k ) - q%(u k ) - L ■ a n ) = 

fc^oo 
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is also true. Now 

lim (vi\u k ) - V2(uj k ); ctf \u k ) - q%' (u k ) - L ■ a n >) = 0. 

Take now convergent sequences of points x k G <J>i n Di,n-\ C H D<i, n -\ 
such that miy(xk) — > (fc — ► oo). Then, according to Lemma 4.26, the projec- 
tions x' of the points x to the (TV— l)-ball phase space of the system necessarily 
fill up a nonempty open piece U' of 5>'„ = 3>j on one hand, which necessarily 
also belongs to fy' n „ = on the other hand. Consequently, $'» coincides 
locally with ^' n /r, a contradiction. The conclusion is that if (II*) holds for the 
(N — l)-ball system, then it does also for the TV-ball system. This completes 
the proof of the inductive step (N — 1) — > N for statement (II*). 

As mentioned before, the base of the induction on N is N = 2 or 3 or 4 
depending on the form of (£, A), and we treat two cases separately. We remind 
the reader of our notational convention to use densely packed indices. 

Case 4/a/l. o\ / a n . In this case the base of the induction is N = 3 or 4. 

The globalization argument provides that the indirect assumption that $i 
and ^ n coincide locally leads to the global containment <J>i C ^> n on Z?2,n-i- 
Then this ought to hold under the side conditions 777,3 = 1714 = 0, too. But this 
is impossible for the following reason: In the case 777,3 = m 4 = the motion 
of the particles {3, 4} has absolutely no effect on the time evolution of the 
{1,2} subsystem. However, the equation defining $1 is obviously sensitive to 
the kinetic data (position and velocity) of the third ball, while the equation 
5(u>) = is written fully in terms of the subsystem {1,2}. The obtained 
contradiction proves the base case of the induction for (II*) in the case u\ 7^ a n . 

Case 4/a/2. o~\ = a n . Now the base of the induction is N = 2. 

Similarly to the previous case, the globalization argument implies the 
containment $1 C * n . This is an assertion completely in terms of the now 
independently evolving two-particle system {1,2}. We will show that in such a 
complexified, two-particle algebraic system the tangential degeneracy S(uj) = 
of the collision a n cannot completely contain the degeneracy \\vf — w^ll 2 = 0. 
(Now the a's only denote the collisions of the two-ball system {1, 2}.) The proof 
will be a simple induction on the positive integer 77. For 77 = 2 the statement 
is obviously true. Assume now that 77 > 2, and the assertion has been proven 
for all smaller values of 77. According to our inductive hypothesis, for a typical 
element u of the a ra -degeneracy 5(lo) = we have that \\vf (u) — v| (<^) || / 0. 
On the other hand, the collision as is known to preserve the quantity \\vi — 1^2 1 1 2 
(the velocity reflection is complex orthogonal), and, therefore, we get that 
||vf(u;) — 7j|(u;)|| 2 / 0, as well. 

The obtained contradiction finishes the discussion of 4/a. 
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4/b. Proof of (III*) by using the inductive hypothesis. Assume now, 
besides the inductive hypothesis, the precondition of (III*). 

The proof of the inductive step (N — 1) — > N is absolutely analogous to 
that given in part 4/a: one has just to substitute $1 with ^i. The verification 
of the basis of the induction will be split into three parts depending on the 
form of (£, .A). In all cases we again assume indirectly that \I>i coincides locally 
with in -D2,n-i- By the globalization/substitution argument from the proof 
of (II*) we then obtain $1 C *S> n to hold globally in I^2,n-i for the base cases, 
too. What we show is that this inclusion leads to a contradiction in all cases. 

Case 4/b/l. o\ ^ a n . In this case the base of the induction is N = 3 or 
4, and the impossibility of the containment *$>i C ^ n for N = 3 or 4 can be 
obtained in the same way as in 4/a/l. 

Case 4/b/2. o\ = o n and there is a k, 1 < k < n, such that cr^ n a n ^ 
and o~k 7^ On- 
Now the base case will be N = 3 and, as explained above, we can assume 
the global containment C ^> n on Z^.n-i- By re-labeling the particles, we 
can assume that a± = a n = {1, 2} and = {2, 3}. 

The global containment C ^> n will naturally be preserved under the 
substitution L = 0, too. This substitution raises the same questions as the 
substitution mjv = 0, and they can be handled analogously. For instance, it 
is easily possible that though (X,«4, r) is minimal from both sides, the scheme 
(£', A', t 1 ) derived from it via the substitution L = is not minimal. Our way 
of taking care of this possibility in the case of the substitution = can 
also be repeated here and, therefore, by our combined indirect and inductive 
arguments, we can again assume that the containment C *$> n now holds 
even in the 3-particle system {1,2,3} (having a time evolution with Property 
(A)). On the other hand, the fact that this 3-particle system {1, 2, 3} also has 
the property L = amounts to a; = for / = 1, . . . ,n. Obviously, we can 
assume that k = 2, i.e., 02 = {2,3}. Case 4/b/2 will be taken care of by the 
following: 

Lemma 4.27. Suppose that the triple (S = (do, . . . , cr n ), 0, r) enjoys Prop- 
erty (A), N = 3, do = o n = {1,2} and, say, o~\ = {2,3}. Then cannot be 
(globally) a subset of ^> n and, therefore, the tangential degeneracies at uq and 
a n (i.e. ^> n ) cannot even locally coincide in Di >n _i. 

Proof. In order to simplify the notation we assume that r = 1/2, and 
o~ n -i / {1,2}. Denote by S* = (a*,..., a*) the symbolic sequence that can 
be obtained from our original sequence (do, . . . , a n ) by discarding all symbols 
{1,2}. The truncated sequence S* consists of a pure block of n(l) {2,3}- 
collisions first, then a pure block of n(2) {1, 3}-collisions etc., up to the closing 
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pure block of n(k) collisions of type {(3 + (— l) fc+1 ) /2, 3}. Consider first 
the following, essentially two-dimensional and completely real, limiting set- 
up: mi = 777-2 = 1, m 3 = 0, ql = 0, q% = e i: q% = -e 2 , = v% = 0, 
v o = (i +£ 2)-i/2.^ ei+£e2 ^ whereei = (1,0,... ,0) G M 1 ', e 2 = (0,1,0,... ,0) G 
IR 1 ', and e > 0, £ — > is a positive parameter of perturbation. By using 
Property (A) of (S, 0,r) and the entire inductive hypothesis of 4.21, we can 
now construct a (S, 0, r) orbit segment 

w = {qf, v?|i = 1,2,3; a = 0,...,n} 

by extending the fixed initial values mi = 777,2 = 1, mz = 0, = 0, q% = ei, 
q\ = — e 2 , v® = v® = 0, 7J3 = (1 + e 2 )" 1 / 2 • (ei + ee 2 ) in such a way that 

(i) At every collision aj = {2,3} (0 < j < n) the smaller root tj is always 
selected out of the two distinct, real roots; 

(ii) At every collision <jj = {1, 3} (0 < j < n) the greater root Tj is always 
selected out of the two distinct, real roots. 

We note here that, obviously, the above constructed orbit segment to is 
a highly degenerate one because of the degenerate {1, 2}-collisions. However, 
we will transform u> into a nondegenerate orbit segment (except one tangential 
collision at <7q) by the upcoming, final perturbation. 

For j = 1, . . . , k we denote by the total number of collisions (i.e. 
together with the {1, 2}-collisions) from <7i up to the closing collision o* h (h = 
n(l) + • • • + n(j)) of the j th X* block. We note that, due to the selection of 
the initial kinetic data, the {1, 2}-collisions have no effect on the motion of the 
third disk. 

An easy geometric calculation shows that the asymptotic formulae 
qf\ £ ) = ei-e 2 - (2n(l) - 1) • (26) 1 / 2 ■ ei + o (e 1/2 ) , 
v f\e) = ei- 2n(l)(2e) 1 / 2 . &2 + o (e 1 ' 2 ^ 
hold. By similar arguments we obtain the asymptotics 



^( 2 ), 



(e) = -e 2 + (2n(2) - 1) • [477(1) • (26) 1 / 2 ] 1/2 . e± + ( £ V 4 ) , 



1 /9 

vl {2 \e) = ei + 2n(2) • [4n(l) • (2e) 1/2 ] • e 2 + (e 1 '^ 

Then a simple induction yields that, in general, the following asymptotic for- 
mulae hold for j = 0, 1, . . . , k: 

(4.28) 

% U \s) = 1 + ( ~ 1)J+1 • ei - e 2 + (-l)^ 2 " ■ ei + o (e 2 "') , 
„^)( £ ) = ei + ( 1)V,~ 2 • e 2 + o (e 2 "') , 
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where the positive constants bj and Cj can be computed by the recursion 60 = 0, 
c = 1, bj = (2n(j) - 1) y/lcj-x, Cj = 2n(j) y/2c~[, j = 1, . . . , k. 

Let us now focus our attention on the almost vertical velocity changes 

AVj(e)=v^\s)-v^- 1 \e) 

and their horizontal components (AVj(e); ei), j = l,...,k. The asymptotic 
formulae (4.28) for the velocities immediately yield the expressions 

(4.29) AVj(e) = (-l) j c j£ 2 ~ j -e 2 + o (e 2 ~ J ) , 

{AV J (ey i e 1 ) = - 1 - C y 1 -> + o(e^), 

j = 1, . . . , k. Here we used the fact that the norm ||f|(e)|| is always 1. 

From the asymptotic formulae (4.29) we conclude that the order of mag- 
nitude of the horizontal velocity change is the highest one in the last pure 
block of collisions B k = (cr*_ n(fc)+1 , . . . , cr*) (actually, s - n(k) = Y?jZ{ n (j))- 
These collisions take place between the third disk and the disk with the label 
i(k) = (3 + (— l) fc+1 ) /2. However, after a perturbation (to be carried out a 
bit later) 1713 ~ 0, G M + (i = 1,2,3), only a certain part of the horizon- 
tal momentum given by the third disk to the i(/c)-disk (during the collisions 
of will be retained by the disk i(k): the rest will be transported to the 
other heavy disk with label 3 — i{k) by the almost tangential collisions {1,2} 
(taking place in between the collisions of B^) with almost horizontal normals 
of impact. This justifies the following definition: 

H : = {q G Z| 1 < q < n(k) and the number of the 

{1, 2}-collisions between o~*_ n ^ +q and a* is even}, 
H:={l,2,...,n(k)}\H. 

By using the method of the proof of the second line of (4.29) we see that 
the order of magnitude of the horizontal velocity change at the q th collision of 
B k {l<q< n{k)) is 

A(2q - l)c k ^- k +o(e 2l - k ). 

Note that the cumulative horizontal velocity change caused by the collisions 

* * 

°s-ra(fc)+l' • • • ' s-n(k)+q 

has the asymptotic formula 4q 2 c k -i£ 21 k ; note how one can deduce the formula 
c k = 2n(k) y /2c k - 1 . 

There are now two possibilities: 
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Case A. £, eH (2(z - 1) + \ [n(k)] 2 . (Note that £, 6Hul? (2<z- 1) = [n(fc)] 2 -) 
In this case the small perturbation 7773(e), 

(4.30) < 7773(e) << e, is very small, depending on e, 

v Q(e) = -v%(e) = 77(e) • e 2 , < 77(e) = r)(m 3 (e)) « m 3 (e) 
and 77(e) is very small compared to 771,3 (e) 

(still mi = "1-2 = 1) serves our purpose: the horizontal component of the 
relative velocity right before the collision a n = {1,2} will be of the order of 
magnitude 



(4.31) 4777 3 Cfc_l 



q£ H q eH 



nl-k 

e 2 + 0. 



We note that the asymptotic formulae (4.28-4.29) will not be destroyed as long 
as the quantities 7773(e) and 77(e) << 7773(e) are small enough depending on e. 
Also note that the collision a n (as well as all other collisions {1,2} amongst 
ao,. . . ,a n ) will be an asymptotically tangential collision with asymptotically 
horizontal normal of impact. These facts immediately follow from the asymp- 
totic formulae (4.29). Hence, in Case A we managed to perturb out of the 
degeneracy \& n from within ^o- 

cwb. £, e ff(2<z-i) = £, e3? (2 9 -i). 

In this case the equation 7771 = 7772 does not suffice: we need to use 
the possibility of selecting different "heavy" masses. First we select and fix 
the positive values mi 7^ 7772, then we carry out the perturbation (4.30) and 
compute the orders of magnitude of the horizontal components Vf 1 , VJ? of 
the velocities of the first and second disks measured right after the collision 
Cn-i = a*s- By applying the collision laws for horizontal "head-on" collisions 
with masses m\ and 7772 we easily obtain the asymptotics 

(4.32) 



2?77i?773 7773 7773 



V? = 2c fc _! [n(k)f 

L7772(777l+ 7772j 777i "72 
Trh a r /,\i2 77717773 nl-fc / ol-k\ 

V 2 h = 4c fc _i 77, (k — — l —l — - • e 2 + e 2 . 

mo 7771 + 7779 V / 



e 2l "+ofe 2 



771-2(7771 + 7772) 

These asymptotics are different, thanks to the inequality mi 7^ m2- Therefore, 
we again managed to perturb out of the degeneracy from within ^0 by 
using the fact that all collisions {1,2} amongst uo, . . . , o n are asymptotically 
tangential collisions with asymptotically horizontal normals of impact. 

Finally, by using the method of obtaining different asymptotics (4.31)- 
(4.32) for the horizontal components of v^ 1 and v^" 1 , we can also achieve 
different asymptotics for the horizontal components of v\ and right before 
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every collision Ui = {1, 2} (i > 0) by appropriately selecting the perturbation. 
In that way all collisions a±, . . . ,a n become nonsingular, and all kinetic data 
of the constructed orbit segment will be real. 

This finishes the proof of Lemma 4.27. □ 

Our last outstanding task in proving (I)-(VII) is the handling of 

Case 4/b/3. o\ = a n (= {1, 2}), and for every k G {1, ... ,n} it is true 
that either au = o~n or a^H a n = 0. 

In this case not all adjustment vectors at with Ok = o n {p < k < n) are 
the same. 

The base case now is N = 2, and according to the indirect assumption 
we have a global containment of the tangential degeneracies: C *$> n in the 
2-ball system {1,2}. The discussion of Case 4/b/3 will be finished as soon as 
we prove the following: 

Lemma 4.33. Suppose that N = 2, the triple (T,,A,t) enjoys Property 
(A), and not all adjustment vectors a±,...,a n of A are the same. Then 
and ^ n cannot locally coincide in D2,n-i- 

Proof. 

1°. We will construct a one-parameter family 

uj(s) G Cl (£, A, t) = O, (se/cR), 

(/ is an open and bounded interval of R) of (S, A, f)-orbit segments in such a 
way that this bunch of trajectories will explicitly show that and ^ n are not 
the same. In fact, we will see that and ^ n intersect each other transversally 
at a large collection of points of intersection. 

2°. Throughout the construction of the family s 6 1} we shall follow 

these guiding principles: 

(a) The dimension u in this construction will be two: The transversality 
of *i and ^ n in a two-dimensional section g*(s) G C 2 C C u , v\{s) £ C 2 C C 
clearly implies the transversality in the general set-up v > 2. 

(b) All data <?*(s) G R 2 , vj(s) G R 2 , mi(s) = m 2 (s) = 1, r = 1/2, 
L(s) = L » 1 will be real, the side length L of the 2-torus being a very large, 
fixed positive constant L. 

(c) Due to the translation invariance of the system, we will only deal with 
the relative motion 

q\s) := q{(s) — q^s) G R 2 , 

v \s) : = „t( a ) _ „t( s ) G M 2 , ||^(s)|| 2 = 1, (s G /). 
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(d) In virtue of the condition Ya=i 1 1°« — °i 1 1 2 > 0; by cutting off a suitable 
starting segment of E (if necessary) , we can assume that = a\ ^ a,2 ■ 

(e) The construction of the family uj(s) will be such that for every k = 
2, 3, . . . , n the polynomial (3.8) has two distinct, nonnegative real roots. We 
will always take the greater one as Tfc. The inductive assumption on the ir- 
reducibility of these polynomials (3.8) over the field ~Kk—i (which are, in fact, 
irreducible, unless they have to be reducible, because of the repetitive adjust- 
ment vectors.) Part of the inductive hypothesis of 4.21 and Property (A) of 
(Y±,A,t) together ensure that the time evolution of u>(s) will follow the alge- 
braic pattern encoded in (T,,A, f). 

(0 

q°(s) := (cos(^o + s), sin(0 o + s)) + (1 - s) ■ (- sin(^ + s), cos(0 o + s)) , 
v°(s) := (- sin(0 o + s), cos(0 o + «)) , 

where the ray {q°(0) + tv (0)| t > 0} passes through the point L • 02- 

The above constructed family [q°(s), v°(s)\ \s\ < Eq, s£ M} is a so-called 
convex family of rays, or convex wave front; see, for example, [B(1979)] or 
[W(1986)]. This means that 

(± q °( S );v%s))=0 1 (±q°(sy±v°(s))>0. 
\ds I \ds as I 

The time evolution of such diverging families of trajectories is discussed in the 
literature; see [W(1986)]. It is shown (see Section 3, especially Theorem 3) 
that our two-dimensional billiard with inner reflections at the circles 

C a := {q G M. 2 \d(q, L ■ a) = l} (a G 1?) 

has the so-called convex scattering property; i.e. every such convex family 
of rays (with not too big curvature, see [W(1986)]) focuses after any inner 
reflection at C a not later than when it reaches the midpoint of the chord of C a 
along which it is traveling after the considered reflection. 

3°. Plainly, the trajectory uj(s) intersects the circle Co tangentially in 
negative time t±(s) = s — 1; that is, the convex family of rays {q°(s), v°(s)\ \s\ 
< £0} is the "synchronized" family of trajectories belonging to ^1. We recall 
that r = 1/2 and ai = <E 1? . 

4°. The inner reflections at the circles C a2 , . . . , C an (yet to be constructed) 
will be in good harmony with the principle 2°(e) above. 

5°. If L is large enough, then there is a small, open subinterval I2 C 
(— £0, £0) (even the closure of I2 will be a subset of (—£0, £0)) such that the 
rays |g°(s) + tv°(s) \ t > 0} intersect the circle C a2 in two distinct points for 
s G I2, while these rays are tangent to C a2 for s G dh- For s G I2 we elongate 
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the above ray up to the second point of intersection with the circle C a2 , and 
this way obtain the orbit segment up to the first collision with that circle. 

6°. Let a,2 = 0,3 = • • • = a k (2) (2 < k(2) < ra), where either a k ( 2 )+i a k(2) 
or k(2) = n. Suppose that k(2) is still less than n. By inspection of the 
geometry of the sequence of inner reflections 02, • • • , &k{2) a * tne circle C a2 it 
is easy to see that the outgoing velocities v k<y2 \s) (s G I2) sweep out the entire 
space S 1 of velocity directions. Therefore, there is an even smaller interval 
I 3 , C1(J 3 ) C I 2 , such that the rays {q k ^(s) + tv k(2 \s)\t > 0} intersect the 
circle C afc(2)+1 in two distinct points for s E Is, while these rays are tangent to 
C afe(2)+1 for s G 9/3. (Here q k( - 2 \s) and v k( - 2 \s) denote the position and the 

outgoing velocity at the reflection &k(2) taking place at time Yli=i T i-) 

For s G I3 we again elongate the above ray up to the second point of 
intersection with the circle C afc(2)+1 , and this way obtain the orbit segment up 
to the first collision with that circle. 

7°. Now we keep going with the iterative construction of the orbit seg- 
ments lo(s) G (l(E,A, f) described above, by also constructing the nested 
sequence of open intervals I2 D I3 D • • • D Ii ■= I , where k(l) = n. The two 
orbit segments {o;(s)| s G dl} are now obviously in the transversal intersection 
of and Thus the lemma is proved. □ 

So far we have settled statements (I)-(VII) of the Main Lemma. 

5. Our only remaining task in the proof of the Main Lemma is to verify 
statement (VIII). Similarly to (II) (III) it is again sufficient to demonstrate 
the slightly weaker assertion (VIII*) formulated in point 4. 

Assume now, besides the inductive hypothesis, the precondition of (VIII*). 
We again use double induction in the same way it was explained before 4/a. 
Also, the proof of the inductive step (N — 1) — ► iV is absolutely analogous to 
that given in part 4/a: one has just to substitute *$> n with $ n . The verification 
of the basis of the induction will be split into two parts depending on the form 
of E. In all cases we again assume indirectly that $1 coincides locally with 
By the globalization/substitution argument from the proof of (II*) we then 
obtain $1 C $ n to hold globally for the base cases, too. What we show is that 
this inclusion leads to a contradiction in both cases. 

Case 5/1 . o\ 7^ a n . In this case the base of the induction is N = 3 or 4, and 
the impossibility of the containment $1 C & n for N = 3 or 4 can be obtained 
in the same way as in the proof of (II*). 

Case 5/2. o\ = a n and there is a k, 1 < k < n, such that u k n a n 7^ and 
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Now the base case will be N = 3, and as explained above we can derive 
the global containment $1 C <£ n - By re-labeling the particles, we can assume 
that o\ = a n = {1, 2} and = {2, 3}. 

The global containment $1 C <3? n will naturally be preserved under the 
substitution L = 0, too. The assumed containment $1 C <£ n now implies that 
even in the 3-particle system {1, 2, 3} (having a time evolution with Property 
(A)) one has a global containment of the ^-degeneracies at o\ and a n . The fact 
that this 3-particle system {1,2,3} also has the property L = amounts to 
ai = for / = 1, . . . , n. Obviously, we can assume that k = 2, i.e. 02 = {2, 3}. 
Case 5/2 will be taken care of by the following analogue of Lemma 4.27: 

Lemma 4.34. Suppose that the triple (S = (o"o, . . . , cr n ), 0, r) enjoys Prop- 
erty (A), TV = 3, (To = o n = {1,2} and, say, u\ = {2,3}. Then $0 cannot be 
{globally) a subset of <& n and, therefore, the 3>- degeneracies at o"o and a n {i.e. 
$0 and $ n ) cannot even locally coincide in Z?i jra _i. 

Proof. The construction used in the proof of Lemma 4.27 will be used 
in this construction with 77(e) = 0. Then, because vf{e) = v®{e) = 0, the 
initial perturbation lies in <J>i. At the same time, in view of (4.31,4.32), for 
the horizontal components V"/ 1 , V 2 h of the velocities of the first and second 
discs measured right after the collision cr n _i = a* we have \\Vi — V 2 h \\ 2 / 
showing that the perturbation necessarily leads out of Thus the global- 
ization/substitution argument provides the claimed transversality of $1 and of 
$ n , the statement of Lemma 4.34. □ 

This finishes the proof of Main Lemma 4.21, as well. □ 

Corollary 4.35. For every triple (S,A,f) with Property (A) and N 
> 3, there exists at least one derived scheme (£', A', t 1 ). 

Corollary 4.36. For combinatorial schemes satisfying Property (A), the 
quadratic polynomial bix 2 + + is reducible over the field K;_i {i.e. Ki_i = 
Kj) only in the case described above in 3.31. 

Proof. Unless we are considering the case described in Remark 3.31, state- 
ment (V) of Main Lemma 4.21 gives the existence of an to G Cl n -i with 5{ui) 
= 0. But then part 3 of the proof of the main lemma provides \K n : K„_i| = 2. 

□ 

Proof of the key lemma. Final induction 

The next lemma establishes the link between the orbit segments of (H,A,t) 
and (£', A', r'). It is a direct consequence of Remark 4.13 (actually, of the fact 
that the reducibility = Kj is inherited by the (£', A', r') -dynamics during 
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the substitution win = 0) and of part (2) of the definition of the derived 
schemes. 

Lemma 4.37. Fix a combinatorial scheme (T,,A,t) satisfying Property 
(A) and an arbitrary scheme (E',A',t') derived from it. Then, for almost 
every choice of the vector of data (q N ,v N ) EC"x C u , 

(1) The subset {(x', q^, v N ,mN = 0)| x' G C h }DD (X, A, f) is open and dense 
in the submanifold |(x' ,q N ,v N ,mN = 0)| x' G C h } (here h = (2u + \)(N — 

i) + i); 

(2) For a nonempty, open set of orbit segments cu' G 0, (X', A', t') there exists 
a trajectory segment to £ Q (X, A, f) for which mjv(u;) = 0, trunc(cj) = to', 
and (q%(id), v° n (lo)) = (q%, v%); 

(3) For an open and dense set of orbit segments u>' £ Q (£', A', t') there exists 
a trajectory segment u> G (l (X, A, f) for which rriNioj) = 0, trunc(cj) = to', 
and (q%(uj), v%(u)) = (q%, v%). 

(Here "trunc" denotes the cancellation of all kinetic data of the -/V th ball with 
zero mass.) 

Proof. Statement (1) is an obvious consequence of Main Lemma 4.21/(1) 
and the Fubini theorem. The second assertion follows immediately from the 
definition of the derived schemes; see also Remark 4.12. To demonstrate (3), we 
select and fix a typical pair (g^, v N ) according to (1) of Lemma 4.37. Intersect 
the open and dense subset of C h - featured in (1) - with the open and dense 
domain D (X', A', r') and obtain the open and dense intersection D' C C h . 
Suppose now that to' £ ft (X', A', t'), and the vector x' = x(lo') of the initial 
variables of to' belongs to D' . (This is an open and dense condition on to' '.) 
Launch the (X, A, r)-dynamics with the initial data (x' , q^, v N , rajy = 0) € 
D(E,A, t) and try to continue the (X, A, f)-dynamics by also preserving the 
time evolution of the given u' E CI (X', A', r'). Remark 4.13 (more precisely, the 
inheritance of the reducibility proved there) and the definition of the derived 
schemes (especially (2) thereof) ensure that the continuation of the (T,,A,f)- 
dynamics is, indeed, possible by the preservation of the given time evolution 
of to'. Hence Lemma 4.37 follows. □ 

The next lemma will use the following: 

Definition 4.38. Suppose that two indices 1 < p < q < n and two labels 
of balls i,j G {1,...,7V} are given with the additional requirement that if 
i = j, then i G U^p+i ^- Following the proof of Lemma 4.2, denote by 
Qi(x),Q2(x), . . . ,Qu(x) (x G c( 2u + 1 ) N + 1 ^) the polynomials with the property 
that for every vector of initial data x G D (X, A, f) and for every k, k = 1, . . . , v, 
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the following equivalence holds true: 

(^3lo £ 0, such that x(u) = x and (v p i (uj)) k = (^^(uj)^ ^ 

Q k (x) = 0. 

Lemma 4.39. Assume that the combinatorial- algebraic scheme (H,A,r) 
has Property (A), see Definition 3.31, and use the assumptions and notation 
of the above definition. 

At least one of the polynomials Q\, . . . , Q u is nonzero, i.e. vf(u) / vj 1 (oj) 
for almost every uj £ Cl (£,«4, r), see also the Corollary 4.7. 

Proof. Induction on the number N > 2. 

1. Base of the induction, N = 2: First of all, by performing the sub- 
stitution L = 0, we can annihilate all adjustment vectors; see the proofs of 
parts 4/a,4/b in Main Lemma 4.21, in accordance with statement (VI) of that 
lemma. Then an elementary inspection shows that for any selection of pos- 
itive masses (mi, w^), indeed, vf(to) ^ almost surely in the section 
Cl (£, A, t, m) of Cl(T,,A,f) corresponding to the selected masses, since any 
trajectory segment of a two-particle system with positive masses and A = 
has a very nice, totally real (and essentially two-dimensional) representation 
in the relative coordinates of the particles: the consecutive, elastic bounces of 
a point particle moving uniformly inside a ball of radius 2r. Therefore, the 
statement of the lemma is true for N = 2. 

We note here that there are merely four points in the entire proof of Key 
Lemma 4.1 where we need to use Property (A) explicitly: The case N = 2, 
i = j in the proof of this lemma, the Lemmas 4.27, 4.33 and 4.34. Whenever 
a collision of type {1,2} is repeated an even number of times with the same 
adjustment vector and zero time slots, the resulting velocity change is iden- 
tically zero, because an even power of the same velocity reflection is applied. 
(See Example I at the end of this section.) 

2. Assume that N > 3 and the lemma has been proved for all smaller 
numbers of balls. By re-labeling the particles, if necessary, we can obtain 

(i) N ^i,N ^ j and 

(ii) if i = j, then the ball i has at least one collision between a p and o q with 
a particle different from N. 

For the fixed combinatorial scheme (T,,A, r), select a scheme {Tj',A',t') 
derived from it. By applying the inductive hypothesis to (£', A 1 , t'), we get 
that for almost every oj' G £l (Z',A', t') the inequality wJV) + v q ~\uj') holds. 
Moreover, by virtue of Lemma 4.37/(1), for almost every (q%, v^) £ C x C 
one has (x' , q° N , v° N , win = 0) G D (£, A, f) for a dense, open subset of x' G C h . 
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Let us fix one such pair (q^, v%). By Lemma 4.37/(2), for a nonempty, open 
set of J G S7(£',„4 / ,t / ) it is true that (x(u/), q%, v%, m N = 0) G D(H,A,t), 
v?{oj') / i>J _1 (u/), and there exists an oj G fi(E,.A, r) with mjv(oj) = and 
trunc(w) = u/. 

The assertion of the lemma now follows for the iV-ball systems since the 
bare existence of a single orbit segment with the above nonequality implies 
that not all polynomials Qi(x) are zero. □ 

Lemma 4.40. Assume that the combinatorial scheme (Y±,A,t) satisfying 
Property (A) is given for the N (N > 3) ball system {1,2, . . . , N} and, fur- 
thermore, there are two integers 1 < p < q < n such that N £ a p C\ a q , N ^ aj 
for p < j < q, and if o~ p = o~ q then we require additionally that there exists an 
index j, p < j < q, such that a p n Oj / 0. Denote by (T,',A', t') an arbitrary 
derived scheme of (S, A, f). 

Denote by P\{x), . . . , P s (x) the canonically determined complex polynomi- 
als, the simultaneous vanishing of which is the consequence of nonsufficiency 
of the (E, A, f) orbit segments in the sense of Lemma 4.2. // all these polyno- 
mials are zero, then all the analogous polynomials associated with the derived 
scheme {T 1 ',A',t i ) are also zero. 

In other words, if all (£, A, f)- orbit segments are nonsufficient, then the 
same is true for all (£' ,A ,r')-orbit segments; see also the corollary at the end 
of the proof of Lemma 4.2. 

Proof. Assume that the assertion is false for the (£', A 1 , r')-orbit segments, 
that is, there exists at least one sufficient (£', A', r^-orbit segment or, equiva- 
lently, by Lemma 4.8, almost all (£', A', f )-orbit segments are sufficient. 

By putting together Lemmas 4.37/(1,2), 4.39 and 4.8, we obtain a suffi- 
cient (S, A, f)-orbit. Indeed, for typical elements x' G D (£', A', r') we have 
that 

(4.41) for all J G Q (£',.4 / ,r') with x{J) = x' 

u)' is sufficient and i% / v i q ' 1 { UJ ')- 

By Lemma 4.37/(1) for almost every {q%,v%) £ C x C it is true that for a 
dense open subset of x' £ C h 

(4.42) (x', q%, v%, m N = 0) G D(S,i,f). 

Thus, for almost every (q%,v^) G C u x C u we have the following: For a 
dense, open subset of x' G C h (4.41) and (4.42) hold. Further, the velocity 
v q N 1 (oj) = v p n (uj) can be chosen in such a typical manner that it is not contained 
by the complex line 

^(u/) + a.(^;V) -<>')) (A G C), 
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so that condition (2) of Lemma 4.8 is satisfied. (The possibility of the desired 
choice for v p n (uj) = v q j^ l (oj) is easily seen if we represent our variables at the 
p th collison; cf. the proof of Lemma 3.19.) 

By Lemma 4.37/ (1, 2), all this implies that for a typical {q%,v^) G CxC 
there exists a nonempty, open set G of orbit segments u/ G £1 (£', A', f'), with 
(x(u/), q^, v^, tun = 0) G D (£, A, r), such that for every oj' G G there exists 
an a; G f2(£,.A,T) satisfying mjv(^) = 0, trunc(oj) = u/, (^(w), f^(oj)) = 
(q%, v%), and, according to Lemma 4.8, uj is sufficient. Hence Lemma 4.40 
follows. □ 

The next, purely combinatorial lemma is the last ingredient of the induc- 
tive proof of the key lemma. 

Lemma 4.43. Define the sequence of positive numbers C(N) recursively 
by taking C(2) = 1 and C(N) = f • max {C(N - 1); 3} for N > 3. Let 
N > 3, and suppose that the symbolic collision sequence E = (<ti, . . . , a n ) for 
N particles is C(N)-rich. Then there is a particle, say the one with label N, 
and two indices 1 < p < q < n such that 

(i) N G cr p n a q , 

(ii) ^US+i^, 

(iii) dp = a q => 3j p < j < q and a p n <7j / 0, and 

(iv) £' is C(iV — \)-rich on the vertex set {1, . . . , iV — 1}. 

ifere, just as in the case of derived schemes, S' is the symbolic sequence 
obtained from S by discarding all edges containing N . 

Proof. The hypothesis on £ implies that there exist subsequences 
Si, . . . , S r of E with the following properties: 

(1) For 1 < i < j ; < r every collision of Ej precedes every collision of Ej. 

(2) The graph of Ej (1 < i < r) is a tree (a connected graph without loop) 
on the vertex set {1, . . . , iV}, and 

(3) r > C{N). 

Since every tree contains at least two vertices with degree one and C(N) = 
f • max {C(N - 1);3} = f ■ b (b is shorthand for max {C(N - 1);3}), there 
is a vertex, say the one labelled by N, such that N is a degree-one vertex of 
Ejm, . . . ,£j( t ), where 1 < i(l) < • • • < i(t) < r and t > b. 

Now t > C(N — 1) and, therefore, (iv) obviously holds. 

Let Opt be the edge of E^) that contains N and, similarly, let o q < be the 
edge of £j( t ) containing the vertex N. Then the fact that t > 3 ensures that 
the following properties hold: 

(i)' N G oy n cy, 

(iii)' a v i = cj g ' =^ 3j, p' < j < q' and cr p ' H / 0, cij / a p /. 
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Let a p , o q (1 < p < q < n) be a pair of edges ay, a q > (1 < p' < q' < 
n) fulfilling (i)' and (iii)' and having the minimum possible value of q' — p'. 
Elementary inspection shows that then (ii) must also hold for a p , a q . Lemma 
4.43 is now proved. □ 

We are now able to prove the key lemma by induction on the number 
N > 2. Indeed, for N = 2 there are no nonsufficient trajectories with initial 
data in D (T,,A,f) and, therefore, the assertion of the key lemma is obviously 
true. 

The inductive step (N — 1) — > N can be obtained by simply putting 
together Lemmas 4.2, 4.8, 4.40 and 4.43. The proof of the key lemma is now 
complete in the complex case. □ 

Finally, since a nonzero polynomial Pi{x) (1 < i < s) takes nonzero values 
almost everywhere on the real space we immediately obtain the 

validity of the real version of Key Lemma 4.1. □ 

We conclude this paragraph by presenting to the reader two interesting 
examples. The first of them shows the absolute necessity of imposing Property 
(A) on the discrete structure (£, A, r), see Definition 3.31, while the second one 
sheds light on the limitations of our algebraic approach developed in Sections 
3-4. 

Example I. Suppose that the collision graph of £ is a tree (a connected 
graph without loop) with the possiblity of repeated edges. (This allows an ar- 
bitrarily large number of consecutive, connected collision subgraphs.) Assume 
further that the combinatorial- algebraic scheme (H,A,t) is given in such a 
way that it has the following property: There exists a sequence 1 = k\ < &2 < 
■ • • < k r < k r+ i = n + 1 such that 

(1) kj + \ — kj is an even number for j = 1, . . . , r, and 

(2) cr p = a kj , a p = a kj , and t p = for kj < p < k j+1 , j = 1,... ,r. 

Plainly, the orthogonal velocity reflections corresponding to the collisions 

CfcjjCTfcj+l) • • • , CTk j+1 -l 

are the same (j = l,...,r) and, therefore, their product is the identity op- 

n k 1 

erator. Thus v'(lo) = (lo) for all i = 1,...,N, j = 1, . . . , r + 1, u G 

f2(£,„4, r). It is an easy exercise then to show that these orbit segments uj 
are very far from being sufficient. Actually, dimc{«i, . . . ,a n } = N — 1 and 
dim c AA(cj) = v + N — 1; see also (4.4). 

Example II. Suppose now that the C(iV)-rich symbolic sequence £ only 
contains the collisions of type {i, i + 1} for i = 1, . . . , N — 1; i.e. the collision 
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graph is the simple path of length N with the allowance of repeated edges. 
Consider the set X of orbit segments lj € Q(E,A,t) with the side conditions 
mi = 7713 = 0. (Of course, rri2 / 0.) It is easy to see that the orbit segments 
from X are never sufficient; the vectors 

(5q° 1 ,...,5q N )={X-(v 1 -v° 2 ),0,...,0) 

(0 / A £ C) are nontrivial neutral vectors. The interesting feature of this 
example is that the set X has two codimensions in O (H,A,f), and the codi- 
mension of X decreases to one at the crucial step of the induction when we 
perform the substitution m\ = 0. (In earlier notation mjv = 0.) 



5. Proof of the Main Theorem 

Richness will be understood throughout in the sense of Definition 2.5 by 
choosing C := C(N) as prescribed by Key Lemma 4.1 (see Remark 4.1/b after 
the key lemma). 

According to the results of [V(1979)], [G(1981)] and [B-F-K(1998)], in 
a semi-dispersing billiard system satisfying the nondegeneracy condition men- 
tioned before our Main Theorem, there are no orbits with a finite accumulation 
point of collision moments. Thus, it is sufficient to concentrate our attention 
on the set M° of orbits containing no singular collision, since it is well-known 
and easy to see that phase points with a singular orbit form a countable union 
of proper submanifolds. 

Theorem 5.1. Consider a system of N (> 3) particles on the L-torus T V L 
(y > 2) satisfying r £ Rq. Let P = {Pi,^} be a given, two-class partition of 
the N particles, where, for simplicity, Pi = {1, . . . , n} and P2 = {n + 1, . . . , N} 
(n < N — 1). Then the closed set 

F + = {x eM: S^'^x is partitioned by P} 

has measure zero. 

This statement would, in fact, allow us to fix C(N) (the number of con- 
secutive, connected collision graphs required for sufficiency) arbitrarily large. 

Proof. Our reasoning is reminiscent of the proof of Theorem 5.1 of 
2 [Sim(1992)-I] irrespective of the fact that the notion of center of mass has 
lost its sense now. The two cases min{n, N — n} > 2 and min{n, N — n} = 1 
can be treated similarly, and thus we only consider the first one. 
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Every point x G M can be characterized by the following coordinates in 
an essentially unique way: 

(1) ir Pl (x) = xi G Mi, (2) irp 2 (x) = x 2 G M 2 , 

(3) C 1 (x)= g n+ i-9i€R", (4) € S"" 1 , 

(5) ||/i(x)|| g R + , (6) ^(x) = ^£? =1 m i t;?( a 0€R + , 

\ 777, f^'fccl 

where = — — — . (Non-uniqueness only arises in choosing C\{x) 

as an arbitrary representative of g n+ i — q± G In what follows the six 

coordinates corresponding to the characterization given above will, in general, 
be denoted by bi,. . . ,bg; thus it will always be assumed that b\ G Mi, 62 G 
M 2 ,b 3 G W,b 4 G S u - 1 ,b 5 ,be G R+. The relation /jl(F+) = will certainly 
follow if we show that for almost every such choice of the b^s 

(5-2) Vb u b 2 ,b 3 ,b 5 A( F +( b ii b 2> b 3,h,h)) = 0, 

where 

F+{b 1 ,b 2 ,b 3 ,b 5 ,b e ) = {y A G S^ 1 : (b 1 ,b 2 ,b 3 ,y 4 ,b 5 ,b e ) G 

and Mb! ,6 2 , 63,65,66 denotes the conditional measure of /U under the conditions 
corresponding to fixing the values of b±, b 2 , 63, 65, b§. (This conditional measure 
is equivalent to the Lebesgue measure on its support.) 

The relation x G F + is equivalent to saying that for every pair % G P\ , j G 
P 2 and every i > 

(5.3) g(ql(x)-q t j (x),L-Z' / )>2r, 

where £?(.,.) denotes the euclidean distance. For simplicity, fix i = 1 and 
j = n + 1. 

Now we will consider the subdynamics corresponding to our two-class 
partition (cf. the "Subsystems, decompositions" part of Section 2 in the paper 
[Sim(1992)-I]), and will denote them, for simplicity, by Si and S 2 , respectively 
(their phase spaces are Mi and M2, of course). 

It is worthwile to note here that — according to the construction, given 
in the introduction, of the configuration space of our hard ball system with 
an arbitrary mass vector — we cannot talk about the absolute configuration 
vector of a particle but rather about the relative position qj(x) — qj(x) for any 
pair of particles, only. Also, it makes sense to talk about the time-displacement 
q\(x) — q®(x) for any 1 < i < N, iet and x G M since 

q\{x) - q°(x) = [ v s i {x)ds (mod L ■ Z"). 

Jo 

Our remarks naturally apply to the subdynamics Si and S2 as well. 
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Now, assuming for simplicity 1 < n < N — 1, a standard calculation, which 
is given for instance in the proof of Lemma 5.8 of [Sim(1992)-I], yields that 

(5.4) 

r(5t rat 

ql+i(x) - Qi(x) = / v s n+ i{x 2 )ds - / uf (xi)ds + Ci(x) + It (mod L ■ Z") 
Jo Jo 

where vf(xi) and vf l+1 (x 2 ) denote the time evolution of the corresponding 
velocity vectors under the subdynamics S\ and S 2 , respectively, and / is the 
relative velocity of the "baricenters" of the second and first subsystems. This 
term appears since in Mi and M 2 the moments of the subsytems are scaled 
to be equal to 0; actually, 

(5.5) ~ I= -M 1; M 2li = - Ml M 2 im h 



M 2 M 2 \\h\ 

where 

n N 

.Mi = ^^mj and M. 2 = raj 

i=l i=n+l 

are the masses of the subsystems. Finally, 

a = s/2Ei(x) -Mi\\h(x)\\ 2 , and 



/ M 2 
P=Jl-2E l (x)- J J ^\\I 1 (x)f 



are the corresponding time scalings. 
Our task is to show that the event 



(5.6) 



for all t > 0, q (^JfJ + /(*). L-IJ^j > 2r 



has measure zero for every fixed b\, b 2 , 63, 65, where 

\ = (M 2 )- 1 (-M 1 -M 2 )\\h\\ 

and 

f(t) : R -> R u 

is an arbitrary fixed function such that f(t) = f(t,xi,x 2 ,Ci, ||/i||, £1). Actu- 
ally, by (5.4), 

/(*)= r< + i{x 2 )ds- FvKxJds + dix). 
Jo Jo 

In (5.6), the canonical meaning of zero measure is that in p^y. 

Denote by £ 2 r,L the lattice of balls of radius 2r centered at points of L-Z". 
Our proof of Theorem 5.1 will be based on the following well-known elementary 
lemma. 
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Lemma 5.7. Fix a vector n G S u ~ 1 for which at least one ratio of co- 
ordinates is irrational. Consider arbitrary hyperplanes H perpendicular to n, 
and denote by B R {z) the {y — 1)- dimensional ball of radius R in H centered at 
z G H. Then, for a suitable j(n) > 

(5.8) liminf inf inf ™™P R (z)n C^) 

R^oo H zeH meas(B R (z)) 

Proof. The lemma easily follows from the fact that the set B R {z) modulo 
L ■ Z u is "uniformly dense" in the torus T£; i.e. it is an e(i?)-dense set in T£, 
where e(R) — ► as R — > oo. This statement is, however, an easy consequence 
of the well-known density of the hyperplane H (mod L ■ Z") in the torus 
T V L = W /L ■ Z u and the compactness of Y" L . Indeed, since being an e-dense 
set is a translation invariant property, the positioning of the foot point z is 
irrelevant. For a given positive e, we choose first a finite e/2-dense set F C T", 
and then for every point y £ F we select a radius R y > with the property that 
Br v {z) is closer to y than e/2. The maximum value R[e) of the finitely many 
radii R y (y G F) entering the problem provides now an e-dense set B R t £ \(z). 
Now select and fix e = r. Observe that at least a fixed, positive percentage 
of the ball B R r r \(z) belongs to £2r,L- However, for every large enough R, also 
a fixed, positive percentage of the ball B R (z) can be covered by disjoint balls 
B R ( r )( z i)- Hence Lemma 5.7 follows. □ 

Assume that statement (5.6) is not true; i.e. the measure of the subset K 
of S y_1 described by (5.6) is positive. Select then and fix a Lebesgue density 
point ft of K with the property that at least one ratio of the components of n 
is irrational. Denote by G £ C S^ _1 the ball of radius e around n. By (5.8) we 
can choose Ro so large that for R > Rq 

. f . f meas (B R (z) C\ C 2 t,l) > jjn) 
H zeH meas (B R (z)) ~ 2 

The set XtG £ can be arbitrarily well approximated by a ball of radius Xte = R 
(R is fixed, R > Rq) in the hyperplane orthogonal to ft through the point 
Xtn + fit) if only t is sufficiently large. Consequently, if R > Ro, then by 
choosing t sufficiently large and at the same time putting e = (At) _1 i?, we 
have 

meas {{\tG £ + /(*)) n C 2 r, L ) > l{n) 
meas ({XtG £ + /(<)) ~ 4 

But this inequality contradicts to the fact that n was chosen as a Lebesgue 
density point of the subset K C S 1 ^ 1 . Hence Theorem 5.1 follows. □ 



Return now to the demonstration of the Main Theorem, restricting atten- 
tion to the subset M°. From Theorem 5.1 it follows that almost every point is 
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C(iV)-rich whatever the value of C(N) may be. Let us fix C(N) according to 
the requirement of Key Lemma 4.1. Denote by R C M° the subset of C(N)- 
rich phase points. By the real version of Key Lemma 4.1, apart from a count- 
able union of proper analytic submanifolds of the outer geometric parameters 
(m, L), for every (T,,A, f) satisfying Property (A), almost every nonsingular 
orbit segment with the combinatorial scheme (H,A,t) is sufficient. Since the 
number of combinatorial schemes is countable, and these schemes of actual 
(real) hard ball trajectories necessarily satisfy Property (A), we can conclude 
that, indeed, almost every orbit in R is sufficient. 

At this point we will use the following simple consequence of Poincare 
recurrence and the ergodic theorem. 

Lemma 5.9 (see [S-Ch(1987)]). If x is a nonsingular orbit segment 
and it is sufficient, then, in a suitable open neighbourhood of x, all relevant 
Lyapunov exponents are nonzero almost everywhere. 

Adding Lemma 5.9 to what has been said before, one also obtains that for 
almost every phase point in R, the relevant Lyapunov exponents of the flow 
are nonzero. Thus our Main Theorem is proved. □ 

Proof of the Corollary to the Main Theorem (Section 1). The positiv- 
ity of the ergodic components is a consequence of the Katok-Strelcyn the- 
ory (cf. [K-S(1986)]), where the underlying assumption is that the relevant 
Lyapunov-exponents of the system do not vanish. (We note that, in fact, the 
Katok-Strelcyn theory is a generalization of Pesin's theory (see [P(1977)] for 
hyperbolic systems with singularities, also containing billiards isomorphic to 
hard ball systems). 

Finally, it is worthwhile to note here that, in fact, on any positive ergodic 
component (i.e. on any ergodic component of positive measure) the Bernoulli 
property of the standard billiard flow (M, {S 1 *}, /x) - L follows from its hyper- 
bolicity and weak mixing, as has been recently shown by Chernov-Haskell and 
Ornstein- Weiss; cf. [C-H(1996)] and [0-W(1998)]. For a brief summary of how 
the K-property follows from ergodicity and hyperbolicity, see, for instance, Sec- 
tion 7 of [Sim(1992)-I] , or [K-S(1986)]. □ 

6. Concluding remarks 

1. Had we been able to strengthen the statement of Key Lemma 4.1 from 
"/or almost every to G Q(E,A, t)", i.e. from "for every uj G (l(T,,A, r) apart 
from a closed algebraic subset of codimension one", to a similar statement 
with an at least two-codimensional algebraic subset, we would have obtained 
the global ergodicity of the system. A crucial step in getting to codimension 
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two in the Key Lemma would be the strengthening of the last claim of Lemma 
4.39 to stating vf (u>) / t>J _1 (w) apart from a subset of codimension two. Both 
lemmas would be implied by the following analogous algebraic properties: All 
common divisors of the polynomials 

(1) Pi,...,P s (in case of Key Lemma 4.1); 

(2) Qi, . . . , Q v (in case of Lemma 4.39) 

are purely homogeneous polynomials of the masses (example II at the end of 
Section 4 illustrates the essential difficulty arising in the inductive proof of 
this claim). In this way one can formulate a hypothetical statement about the 
Boltzmann- Sinai ergodic hypothesis: If this stronger form of Key Lemma 4.1 
is true, then, in the Main Theorem, we can also claim the K-property instead 
of just full hyperbolicity of the system. 

2. The system oi N = 2 hard discs (i.e. v = 2) with arbitrary masses 
mi, rri2 was studied by Simanyi and Wojtkowski [S-W(1989)]. They considered 
the same dynamics as here, but in the larger phase space Q x S 1 . By using the 
fact that this system is an isometric T 2 -extension of the standard billiard ball 
flow (M, S*, /^)( mii m 2 ), which is a planar dispersing billiard, and was known 
to possess the Bernoulli property since [G-0(1974)], they showed that this 
extension is also a Bernoulli flow as long as the ratio m\/m2 is irrational. 
If one thoroughly studies the steps of that proof, he realizes that, once the 
difficulties outlined in the above remark have been resolved, the methods of 
the present paper and [S-W(1989)] together with the results of [C-H(1996)] 
and [0-W(1998)] prove even the Bernoulli property of the flow on the extended 
phase space Q x£ (without the factorization see the introduction), provided 
that not all ratios mi/rrij are rational. (The latter condition means precisely 
that the "center of mass" is defined modulo a dense subgroup of T".) 

3. In [Sz(1994)], the second author obtained necessary and sufficient con- 
ditions for the K-property of the so-called "orthogonal cylindric billiards." 
The subclass of cylindric billiards, introduced in [Sz(1993)], within the family 
of semi-dispersing billiards, is of particular interest since it is this class for 
which the formulation of general — and at the same time in some sense con- 
structively verifiable — necessary and sufficient conditions of the K-property 
seems possible. (We have seen in the proof of Lemma 2.1 that the system of N 
balls with an arbitrary mass vector is also isomorphic to a cylindric billiard.) 
Indeed, such a condition is found and shown to be necessary for the systems of 
N hard balls of identical masses in a forthcoming paper of the present authors 
[S-Sz(1997)]. One apparent difficulty in establishing the sufficiency of this con- 
dition is the fact that so far no characterization of the neutral subspace of a 
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trajectory segment of a general cylindric billiard is known (in the case of hard 
ball systems such a characterization was obtained by using the CPF, whose 
derivation, on the other hand, used the conservation of momentum, a property 
present in hard ball systems but absent in general cylindric billiards). 
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